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ABSTRACT 

We present the results of new near-IR spectroscopic observations of passive galaxies at z ^ L4 in a concen- 
tration of 5z^-selected galaxies in the COSMOS field. The observations have been conducted with Sub- 
aru/MOIRCS, and have resulted in absorption lines and/or continuum detection for 18 out of 34 objects. This 
allows us to measure spectroscopic redshifts for a sample that is almost complete to /^ab = 2L COSMOS 
photometric redshifts are found in fair agreement overall with the spectroscopic redshifts, with a standard de- 
viation of ~ 0.05; however, ~ 30% of objects have photometric redshifts systematically underestimated by 
up to ~ 25%. We show that these systematic offsets in photometric redshifts can be removed by using these 
objects as a training set. All galaxies fall in four distinct redshift spikes at z = 1.43, 1.53, 1.67 and 1.82, 
with this latter one including 7 galaxies. SED fits to broad-band fluxes indicate stellar masses in the range 
of 4-40 X IO^^Mq and that star formation was quenched ~ 1 Gyr before the cosmic epoch at which they 
are observed. The spectra of several individual galaxies have allowed us to measure their HJp indices and the 
strengths of the 4000 A break, which confirms their identification as passive galaxies, as does a composite 
spectrum resulting from the coaddition of 17 individual spectra. The effective radii of the galaxies have been 
measured on the COSMOS HST/ ACS /psMW-band image, confirming the coexistence at these redshifts of pas- 
sive galaxies which are substantially more compact than their local counterparts with others that follow the 
local effective radius- stellar mass relation. For the galaxy with best S/N spectrum we were able to measure a 
velocity dispersion of 270 ±105 km s~^ (error bar including systematic errors), indicating that this galaxy lies 
closely on the virial relation given its stellar mass and effective radius. 
Keywords: galaxies: evolution — galaxies: formation — galaxies: high-redshift 



1. INTRODUCTION 

Galaxy formation and evolution in the cosmological con- 
text is one of the leading themes in observational cosmology 
today, with much effort being currently dedicated using all fa- 
cilities both on the ground and in space. The goal is to map 
the evolution of the population of galaxies as a function of 
redshift (cosmic epoch) and environment, with each popula- 
tion being described by distribution functions of mass, star 
formation rate (SFR), morphology, AGN activity, etc. This is 
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^ Based on data collected at the Subaru telescope, which is operated by 
the National Astronomical Observatory of Japan. (Proposal IDs: S09A- 
043andS10A-058). 



accomplished through a full multiwavelength approach, cov- 
ering the whole electromagnetic spectrum, so to ensure the 
most robust characterization of each individual galaxy. This 
growing body of observational evidence has started to provide 
enough data for a fully empirical reconstruction of the main 
trends in the evolution of galaxies, such as their mass growth 
via star formation and merging, and the quenching of their 
star formation activity, t urning the m into passively-evolving, 
early type galaxies (e.g. JPeng et al . 2010b, 201 1). In parallel, 
these observational constraints are continuously fed into the- 
oretical simulations, forcing them to match closer and closer 
the picture that is progressively emerging from observations. 

Quenching of star formation is perhaps the most dramatic 
event in the evolution of galaxies, turning them into passively- 
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evolving galaxies (PEGs). Hence, properly mapping this 
transformation across cosmic times and in a range of envi- 
ronments must be central in our attempt to understand galaxy 
evolution. Indeed, massive PEGs represent the culmination 
of the galaxy evolution process and are known to be the 
most strongly clustered galaxies loc ally as well as at z ~ 2 
(iKong et al.ll2QQ6l : lMcCracken et al.l l201Q). Besides marking 
the highest density peaks, PEGs (ellipticals) are extremely in- 
teresting on their own. In the local Universe PEGs (includ- 
ing bulges) encompass almost 60% of the total stellar mass 
(^Baldry et al. 2004), and therefore the formation and evolu- 
tion of this type of galaxies are central to the broader problem 
of galaxy evolution in general. 

A most extended and coherent effort to map galaxy 
evolution up to z y 3 is certainly the COSMOS project 
(IScoville et al.ll2007h with full multiwavelength imaging cov- 
erage of an equatorial 2 square degree field. Besides the imag- 
ing campaign extending from the X-rays to radio, a major 
spectroscopic effort over this field has been conducted with 
the Very Large Telescope (VLT), ta rgeting galaxies all the 
way to z - 3 (the zCOSMOS project. ^Mj^^^oMK 120091 
and in preparation). As a result, spectra and redshifts have 
been obtained for ~ 20, 000 galaxies up to z — 1 .2 over the full 
COSMOS field (zCOSMOS-Bright) and for - 5,000 galax- 
ies at 1 .4 < z ^ 3 over its central square degree (zCOSMOS- 
Deep). 

Thus, mapping the evolution of PEGs as a function of 
mass, redshift and environment has been fairly well accom- 
plished up to z ~ 1, especially as part o f the zCOSMOS- 
Bright project (e.g., Scarlata et"aD 120071; iBolzonella et al. 
12010: lPozzettTetal]l2010D . and there is a claim that the most 
massive PEGs are virtually all in place by z ~ 1 (ICimatti et al.l 
120061: iPozzettiet ani2010l) . However, their number density 
rapidly falls at higher redshifts as suggested by the number 
density of PEGs at (z) = 1.7 being only 20% of that of 
local ellipticals of similar m ass (Kong et al. 2006 , see also 
Daddi et al.l l2005l: iBrammer et al. 2009 2011: Tavlor et al. 
20091: IWilhams et al] [20091: lllbert et al.l IMol: iCassata et al. 

However, much of the current statistics of z ^ 1.4 PEGs re- 
lies on photometric redshifts, which are not accurate enough 
to properly map the density field at these high redshifts. This 
is especially the case for the highest density regions, where 
environmental effects are presumably strongest, and likely 
to be preferentially inhabited by PEGs which are the most 
strongly clustered population at 1 .4 < z ^ 3 ( Kong et al. zOOo; 
^cCracken et al. 2010). As a consequence, it is not currently 
possible to measure the fraction of red (quenched) galaxies 
at z ^ 1 . 4 as a func t ion of o verdensity, as done at lower red- 
shifts by lPeng et al.l (1201 Ob|) . In fact, due to instrumental lim- 
itations the high-redshift component of the zCOSMOS sur- 
vey has been designed to observe only star-forming galaxies 
(down to ^AB ^25), and therefore none of the ~ 4, 000 z > 1 .4 
passively evolving galaxies brighter than ^ab ^ 23 was in- 
cluded among the targets. Thus, with zCOSMOS-Deep alone 
one cannot study this important galaxy population, hence one 
cannot properly map the highest peaks of the density field at 
these redshifts. 

Beyond z — 1.4, only a handful of such gal axies have been 
observed spectroscopicallv (^e.g.. ICimatti et al . 2004, 2008; 
McCarthy et~ani2004l: iDaddi et al.ll2005L seelDamjanov et al, 
20111 for a recent compilation), and their stellar populations 
dated. Except for such rare cases, the stellar population prop- 
erties of high-z PEGs have been explored only via their broad- 



band spectral energy di stribution (SEP), or w ith intermediate- 
width passbands (^e.g.. IWhitaker et ani201 1|) . Besides locat- 
ing individual PEGs relative to the density field, full spec- 
troscopic coverage is instead essential to alleviate the age- 
metallicity degeneracy and to get dynamical information, in 
particular in the rest- frame optical where features are exten- 
sively used and calibrated for this purpose in local elliptical 
galaxies. 

At redshifts below ~ 1.4 the strong spectral features that 
allow one to measure the redshift of PEGs are the Ca 11 H&K 
lines and the 4000 A break, but by z = 1.4 these features 
move to the near-IR. The only useful spectral feature for opti- 
cal spectroscopy then remains the complex of rest-frame UV 
absorption lines of neutral and ionized iron and magnesium 
in the 2650 < Arest(A) < 2850 spectral range, where the con- 
tinuum of PEGs is very faint. Nevertheless, this feature has 
been used to derive redshifts and stellar ages for most of the 
few passive galaxies at z > 1.4 with spectroscopic informa- 
tion (see references above). However, reaching the extremely 
faint UV continuum of a passive galaxy requir es extremely 
long integrations, which in the case of Cimatti etaH (l2008h 
ranged from 30 to 60 hours per object. 

Alternatively, one can try to follow the Ca II H&K lines 
and the 4000 A break in the near-IR, where the continuum 
flux is much higher (but also the background). This was at- 
tempt ed with lon g- slit spectroscopy with GNIRS on Gem- 
ini bv iKriek et all (|2006, 2008, 2009) resulting in 17 spec- 
troscopic redshifts of PEGs at z — 2.3 judged from the ab- 
sence of emission lines. Multi-object spectroscopy in the 
near-IR offers decisive advantages over optical or single- 
channel near-IR spectroscopy of high-redshift PEGs, and the 
Multi-Object InfraRed Camera and Spectrograph (MOIRCS; 
llchikawa et al. 2006; Suzuki et al. 2008) mounted on the Sub- 
aru telescope offered an unique opportunity amongst 8-1 Om 
class telescopes. In 2006 we then started a project exploit- 
ing this opportunity and, after a few unsuccessful runs, useful 
data have started to accumulate from 2009 on and the results 
are presented in this paper. 

The paper is organized as follows. In Section |2] we describe 
our dataset; redshift identifications are presented and com- 
pared to photometric redshifts in Section |3] whereas the pres- 
ence of several PEG redshift spikes is discussed in Section SI 
In Section [5] a composite spectrum corresponding to ~ 140 
hours of integration is presented. The properties of the stellar 
populations of the program galaxies are investigated in detail 
in Section[6l In Section[7]the structural and kinematical prop- 
erties of these high-z PEGs are derived and compared to those 
of other high-z PEGs with spectroscopic redshifts and of lo- 
cal early-type galaxies. The number and mass densities of our 
PEGs are discussed in Section [H Finally, we summarize our 
main results and conclusions in Section [9l 

Throughout this paper, we assume a flat cosmology with 
Hq = 70km s"^Mpc"\ = 0.3, and Qa = 0. 7. Photomet- 
ric magnitudes are expressed in the AB system ( Oke & GunnI 
1983). 

2. SAMPLE, OBSERVATIONS, AND DATA REDUCTION 

2.1. Sample Selection 

Our MOIRCS spectroscopic target s are p rimarily BzK- 
selected candidate PEGs (pBzKs, D addi et alJf20 04) extracted 
from the ^-band selec ted catalog for the COSMOS field 
(iMcCracken et al.ll201Cih . Among its ^ 4, 000 pBzKs down to 
K < 23, we selected 34 pBzKs with 19.4 <K <21.9 among 
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the brightest objects that satisfied as many of the following 
additional criteria as possible: 

• no d etection at 24 jim in the COSMOS Spitzer /MIFS 
data ('Sander s et aL[l2QQ7l: iLe Floc'h et ani2009h down 
to 80/iJy, corresponding to a ^ 3cr detection limit, to 
exclude star-forming galaxies and AGNs. About 3,100 
of 4,000 pBzKs satisfy this criterion. At z ^ 1 .4, where 
the pBzKs are expected to be located, this flux limit 
corresponds to a SFR of roughly 50Mq yr~\ which still 
fairly high. As an exception to this criterion, one pBzK 
object with higher 24 /im flux was also included in the 
target list; 

• a higher 4.5 /im flux compared to that at 3.6 /im in the 
Spitzer/IRAC bands, i.e., m^^e-^A.s > 0, which is sat- 
isfied by about 2,900 pBzKs. The IRAC photometry 
is ta ken from the publ icly available S-COSMOS cata- 
log (ISanders et al.ll200 7). A red m3.6 -m4.5 color makes 
likely that the object is at z > 1.4, hence enabling us to 
access major rest- frame optical absorption features.; 



• photometric redshift ffrom lllbert et al.ll2009|) Zphot > 1 .4 
to detect 4000 A break in the observed wavelength 
range. About 2,400 pBzKs have Zphot > 1-4. Sev- 
eral pBzK objects with Zphot slightly below 1.4 were 
also included to account for photometric redshift errors. 
About 30% of our selection have Zphot < 1-4; 

• high projected concentration of pBzK galaxies over 
each MOIRCS field of view (FoV), to make sure that 
the highest possible number of suitable targets is ob- 
served for each telescope pointing, thus exploiting 
the multiplex of the instrument. Bright pBzKs are 
strongly clustered whereas their surface number density 
on the sky is relatively low (r^ Mpc a nd H 0.4 
arcmin"^ to the K = 22 limit, respectively; iKong et all 
120061: iMcCracken et al.' "2010^. Besides maximizing 
the yield, targeting pBzK concentrations allows us to 
check whether such concentrations are mere statistical 
fluctuations in surface density, or real physical struc- 
tures. Thanks to the wide 2 square degree field of COS- 
MOS, several concentrations with at least 15 pBzKs 
with K<22 per MOIRCS field of view were identified; 

• region s including some of the 12 brightest pBzKs stud- 
ied by Mancini et al.l (2010), hence ensuring the bright- 
est possible continuum emission and the detection of 
absorption features with as high a S/N as possible. 

Considering the above constraints, w e selected a FoV con- 
taining the object 254025 at z = 1.82 ("Onodera etal."2010b 
which is one of the 12 brightest galaxies in Mancini et al] 
(t2010), and observed it in 2009 finding a possible overdensity 
of pBzKs with spectroscopic redshift at z = 1.82 (see Section 
|4]). Motivated by this finding, we then selected two more ad- 
jacent FoVs to cover a possible large scale structure around 
it as well as to maximize the number of the brightest pBzKs 
in each mask. Thus, the three selected MOIRCS FoVs cluster 
around a concentration of bright pBzKs, but the selected fields 
lie ou tside the central s quare degree covered by zCOSMOS- 
Deep ("Li lly et al.l[2QQ7h . 

Figure [T] shows the projected density map of pBzKs in the 
COSMOS field with K < 22, along with the position of the 
three MOIRCS pointings for which the data presented here 



have been taken. Besides 254025, two other bright pBzKs 
(namely 217431 and 307881) were also included in the sam- 
ple of 12 ultra-massive PEGs whose struct ural parameters 
have been measured by iMancini et al.l (120101) . Together with 
254025, these are the three brightest objects in our sample, 
and therefore are listed at the top of the tables in this paper. 

One MIPS 24 /im detected object (313880, see Section 
13.31) with /24 = 130/iJy was included because of geometrical 
constraints there is no alternative pBzK object that could be 
placed in the MOIRCS mask. The Hst of the 34 pBzK targets 
along with their photometric redshifts and magnitudes are re- 
ported in Table ll] Note that all the pBzKs down to ^ = 21 
included in the MOIRCS explored area have been included 
in the target list. Figure d displays the BzK plot for the COS- 
MOS field (adapted from McCracken et al.ll2010l) with the ob- 
served pBzK targets shown with red symbols. These same 
objects are also shown in the (z-K)-K color-magnitude dia- 
gram displayed in Figure [3] 

After including in the MOIRCS masks the highest possi- 
ble number of pBzKs, the residual fibres were used to ob- 
serve star-forming EzK-selcctcd galaxies (sBzKs, shown as 
blue symbols in Figures [2] and [3]) . The results for these sBzKs 
will be presented and discussed in a future paper. 
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Figure 1. The density map of the pBzK galaxies in the COSMOS 2 squar e 
degree field covered by the K-bmd image from McCracken et alJ (I201QD . 
Black and white dots show the location of the ~ 2800 pBzKs down to ^ = 22 
and the spectroscopically observed pBzKs, respectively. White rectangles 
and circles show the FoVs covered by MOIRCS observations (the mask IDs 
are B, A, and C from north to south, respectively). Slits were placed in the 
overlapping regions between a circle and the corresponding rectangle. Back- 
ground shows the surface density of pBzKs as calculated from 10-th nearest 
neighbour and normalized by an average surface density of 0.41 arcmin"^. 



2.2. Observations 

The near-IR multi-object spectroscopic observations have 
been carried out with Subaru/MOIRCS. The imaging FoV of 
MOIRCS is 7 X 4 arcmin^ with 2 detectors (channel 1 and 2, 
respectively) and slits can be placed within the intersection 
with the 6 arcmin diameter circular region as illustrated in 
Figure [T] The low resolution zJ500 grism was used with 0.7 
arcsec width slits, which covers about 9000-17500 A with a 
resolution R 500. The slit length ranges from 8 to 14 arcsec, 
adequate to subtract the sky background. 

Six masks in total, two masks per FoV, were observed dur- 
ing six separate observing runs as reported in Table |2] We 
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Table 1 

Photometric properties of observed pBzKs 



ID RA Dec B z K B-z z-K Zphot MasklD^ Exp. Time 

(deg) (deg) (mag) (mag) (mag) (mag) (mag) (min) 



Identified 



254025 


150.6187115 


2.0371363 


25. 15 ±0.03 


22.19 ±0.02 


19.41 ±0.01 


3.23 ±0.04 


2.78 ±0.02 


1.7407 


A 


490 


217431 


15U.oo4o939 


1.9497545 


26.11 ± U.IU 


22.07 ± 0.03 


19.43 ± 0.01 


4.32 ± 0.10 


o z:o 1 f\ r\o 

2.63 ± 0.03 


1.2457 


C 


450 


307881 


150.6484873 


2.1539903 


26.47 ±0.13 


22.06 ±0.03 


19.46 ±0.01 


4.68±0.13 


2.60 ±0.03 


1.4201 


B 


550 


233838 


150.6251048 


1.9889180 


27.44 ±0.29 


22.98 ±0.05 


19.95 ±0.01 


4.73 ±0.29 


3.03 ±0.05 


1.8428 


C 


450 


277491 


150.5833512 


2.0890266 


26.31 ±0.10 


23.08 ±0.06 


20.11 ±0.01 


3.50±0.12 


2.97 ±0.06 


1.7936 


A 


490 


313880 


150.6603061 


2.1681129 


25.97 ±0.07 


22.88 ±0.05 


20.27 ±0.02 


3.36 ±0.09 


2.61 ±0.05 


1.3430 


B 


550 


250093 


150.6053729 


2.0288998 


27.21 ±0.25 


23.57 ±0.09 


20.59 ±0.02 


3.90 ±0.27 


2.98 ±0.09 


1.4234 


A 


490 


263508 


150.5677283 


2.0594318 


26.77 ±0.14 


23.14 ±0.06 


20.61 ±0.02 


3.90±0.15 


2.53 ±0.06 


1.5658 


A 


490 


269286 


150.5718552 


2.0712204 


26.73 ±0.14 


23.14 ±0.06 


20.64 ±0.02 


3.85±0.15 


2.50 ±0.06 


1.6654 


A 


490 


240892 


150.6432950 


2.0073169 


> 27.90 


23.40 ±0.08 


20.69 ±0.02 


>4.67 


2.71 ±0.08 


1.5474 


C 


450 


205612 


150.6542714 


1.9233323 


26.32 ±0.09 


23.33 ±0.06 


20.70 ±0.02 


3.26±0.11 


2.64 ±0.06 


1.2932 


C 


450 


251833 


150.6293675 


2.0336620 


26.36 ±0.12 


23.07 ±0.06 


20.52 ±0.02 


3.56±0.13 


2.55 ±0.06 


1.1874 


A 


490 


228121 


150.5936156 


1.9754018 


27.75 ±0.38 


23.73 ±0.10 


20.73 ±0.02 


4.29 ±0.39 


3.00±0.10 


1.7198 


C 


450 


321998 


150.7093826 


2.1863891 


27.50 ±0.32 


23.63 ±0.09 


20.77 ±0.02 


4.13±0.33 


2.87 ±0.09 


1.4320 


B 


550 


299038 


150.7091894 


2.1369001 


> 27.90 


24.02 ±0.13 


20.97 ±0.02 


>4.15 


3.05±0.13 


1.7667 


B 


550 


209501 


150.6645174 


1.9325604 


> 27.90 


24.17±0.14 


21.56 ±0.04 


>4.00 


2.61 ±0.15 


1.3413 


C 


450 


253431 


150.6408360 


2.0378335 


> 27.90 


24.17±0.13 


21.27 ±0.03 


>4.00 


2.89±0.13 


1.5583 


A 


490 


275414 


150.5822420 


2.0857211 


> 27.90 


23.76 ±0.12 


20.80 ±0.02 


>4.22 


2.96±0.12 


1.4332 


A 


490 



No detection 



222961 


150.6455001 


1.9637953 


> 27.90 


24.76 ±0.25 


21.69 ±0.05 


>3.41 


3.07 ±0.25 


1.4755 


C 


450 


229536 


150.6724126 


1.9780668 


26.58 ±0.14 


24.08 ±0.16 


21.56 ±0.05 


2.75 ±0.21 


2.52±0.17 


1.6066 


C 


450 


243138 


150.5947944 


2.0130046 


27.46 ±0.27 


24.65 ±0.23 


21.93 ±0.06 


3.08 ±0.35 


2.72 ±0.24 


1.9387 


A 


280 


251051 


150.5881329 


2.0320317 


27. 10 ±0.20 


24.21 ±0.14 


21.50 ±0.04 


3.15±0.24 


2.72 ±0.15 


1.4514 


A 


280 


255465 


150.6336861 


2.0423782 


27.35 ±0.24 


24.43 ±0.16 


21.66 ±0.05 


3.19±0.29 


2.77 ±0.17 


1.8805 


A 


280 


258867 


150.5698729 


2.0491841 


> 27.90 


24.77 ±0.25 


21.62 ±0.05 


> 3.28 


3.15±0.25 


2.3654 


A 


280 


260120 


150.6126504 


2.0523735 


> 27.90 


23.82±0.10 


21.31 ±0.03 


>4.35 


2.51 ±0.10 


0.9670 


A 


280 


268884 


150.5525002 


2.0719597 


> 27.90 


24.34 ±0.17 


21.31 ±0.03 


> 3.78 


3.03±0.17 


1.6244 


A 


280 


273534 


150.5994479 


2.0819904 


27.50 ±0.27 


24.52 ±0.20 


21.49 ±0.04 


3.25 ±0.34 


3.03 ±0.20 


1.5801 


A 


280 


281751 


150.5802008 


2.1000786 


27.03 ±0.18 


23.72 ±0.09 


21.10±0.02 


3.58 ±0.20 


2.62 ±0.09 


1.1811 


A 


280 


305677 


150.6343972 


2.1488364 


27.11 ±0.33 


23.04 ±0.09 


20.52 ±0.02 


4.34 ±0.34 


2.52 ±0.09 


1.3697 


B 


340 


315704 


150.6874049 


2.1729448 


> 27.90 


24.10±0.14 


21.37 ±0.04 


>4.05 


2.74 ±0.15 


2.5197 


B 


550 


316338 


150.7264016 


2.1721785 


> 27.90 


24.08 ±0.15 


21.44 ±0.04 


>3.95 


2.64 ±0.16 


1.1895 


B 


340 


321193 


150.6749498 


2.1850247 


27.45 ±0.25 


24.34 ±0.16 


21.54 ±0.04 


3.37 ±0.30 


2.80±0.16 


1.5988 


B 


550 


322048 


150.6514093 


2.1861803 


27.18 ±0.27 


24.03 ±0.14 


20.95 ±0.03 


3.42 ±0.30 


3.08±0.14 


1.5022 


B 


340 


325564 


150.6824812 


2.1934154 


> 27.90 


23.57 ±0.10 


20.72 ±0.02 


>4.31 


2.85±0.10 


1.3220 


B 


550 



Note. — Upper limits are 3cr. 
^ The mask IDs are given in Table |2] 

replaced targets with alternative ones when neither contin- 
uum emission nor emission lines were seen in the first observ- 
ing runs, but the majority of the targets are common between 
the masks targeting the same FoVs (cf. Table [2]). Sequences 
of 600-1200 sec integrations were made in a standard two- 
position "AB" dithering pattern separated by 2 arcsec. At the 
beginning and/or end of the nights AOV-type standard stars 
were observed for flux calibration (i.e., atmospheric absorp- 
tion and instrument response) with the identical instrumental 
setup as the targets and with similar airmass as for the COS- 
MOS field observations. 

In 2009, observations were affected by cloudy weather and 
~ 1.2 arcsec seeing whereas in 2010 we took advantage of 
clear nights and ~ 0.4-1.0 arcsec seeing conditions. Each 
FoV was integrated for 7-9 hours (about half of this if tar- 
gets were replaced between runs). Table |2] summarizes our 
observations. 

2.3. Data Reduction 

The data were red uced with the MCSMDP pipelineEl 
(|Yoshikawa et al.|[2Q10l) and with custom scripts. The data are 
first flat-fielded using dome-flat frames taken with the same 

^ |http : / /www . nao j . org/ Qbserving/DataReduct ion/| 



configuration of the science frames. Bad pixels were then 
removed using the bad-pixel map bundled in MCSMDP and 
cosmic ray hits were removed using the pair of images in the 
dithering. The counts of the rejected pixels were replaced by 
linear interpolations from neighboring pixels along the spa- 
tial direction. After removing bad pixels and cosmic rays, the 
sky background was removed by subtracting the "B" frame 
from the 'A" frame. Then the distortion of the detectors was 
corrected by using the same coefficients used in the MOIRCS 
imaging data reductiorOI. 

Each frame was then rotated to correct the tilt of the spectra 
on the detector by making the edges of each slit almost paral- 
lel to the x-direction (dispersion direction) of the detector. We 
applied a rotation angle of ~ -0.56 and ~ -0.73 degrees for 
the channel 1 and 2, respectively. 

Individual spectra were cut from the frames and wave- 
length calibration was made for each extracted 2-dimensional 
spectrum based on th e location of the OH-airglow lines 
(iRousselot et al.l [200Qh present in the observed wavelength 
range. The uncertainty of the wavelength calibration is typ- 

^ The distortion coefficients are available as an 
IRAF/geotran format bundled with a MOIRCS imag- 
ing data reduction pipeline provided by Ichi Tanaka from 

http : / /www . nao j . org/staf f /ichi/MCSRED/mcsred . html 
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Table 2 

Observation Logs 



Mask ID 


Date 


Exp Time 
(min) 


Number of pBzKs^ 


Number of sBzKs^ 


A 


13, 14 Mar 2009 


280 


16 


16 




29, Mar 2010 


210 


8(8) 


20 (8) 


B 


6 Feb 2010 


340 


10 


17 




1 Apr 2010 


210 


7(7) 


15(7) 


C 


7 Feb 2010 


330 


8 


15 




21 Feb 2010 


120 


8(8) 


14(9) 



Note. — Column 1: ID of the three MOIRCS pointings (masks); Column 2: observing date; 
Column 3: exposure time; Column 4: number of pBzKs in the mask; Column 5: number of sBzKs in 
the mask 

^ The number of common objects between two runs with the same MOIRCS pointing are indicated in 
parentheses. 




12 3 4 

(5-z)ab 



5 6 



Figure 2. The BzK 2-color diagram for objects in the COSMOS field. The 
/sT-selected objects are shown in gray scale as a background and the pBzKs in 
our spectroscopic sample are shown in red. Red squares and circles show 
pBzKs with and without continuum detection, respectively. Blue circles 
represent star-forming BzK-^t\tcit6. galaxies that were also included in the 
MOIRCS masks which will be described in a forthcoming paper. 



ically one half the pixel size, i.e., ~ 2.5 A. Then the frame 
was transformed to align the sky lines along the j-direction 
(spatial direction) and to make wavelength a linear function 
of pixel coordinates. Because sky background residuals are 
expected due to the time variation of sky brightness (in par- 
ticular in the OH-line strength) we performed an additional 
residual sky- subtraction by subtracting at each wavelength a 
mean of the spatial pixels outside those occupied by the ob- 
jects. 

The spectra of the AOV-type stars were processed in the 
same manner as for the targets, and the 2D frames were 
co-added. Then the ID spectra of the standard stars were 
extracted by using the apall task in IRAK The system 
response curves including atmospheric, telescope, and in- 
strumental throughputs were obtained by dividing the ob- 
served spectra of the standard stars by spectra from the stel- 
lar spectral library of Pickles (1998). After the residual sky- 
subtraction, the 2D spectra of the targeted galaxies were flux- 
calibrated by dividing them by the system response curve and 
co-added with appropriate offsets derived either from the lo- 
cation of reference stars or bright compact objects to which 
slits had been assigned or from the fits header information. 



gM/Mo 12.0/ 11.8/ 11.6/ 11.4/ 11.2/ 11.0/ 10.8/ 10.6/ 10.4/ 

/ / / / // // / 

/ / / / / T / /t / 




Figure 3. The iz-K)-K color-magnitude diagram for objects in the COS- 
MOS field. Symbols are the same as in Figure |2| Dashed lines show constant 
stellar mass derived by using Equation (6) and (7) of Daddi et al.. (2004), i.e., 
log(M*/lOiiM0 =-0.4(/rt°t-21.38) + 0.218(z-/s:-2.29), converted to the 
Chabrier IMF adopted in this paper. 

In performing these co-additions we applied weights propor- 
tional to the exposure time and inversely proportional to the 
seeing size and the square of the S/N of each individual 2D 
spectrum. 

Finally, for the objects with continuum detection the ID 
spectra were extracted from the stacked 2D frames by using 
the apall task in IRAK The ID spectra from the different 
observing runs were extracted separately, and then the pairs 
of ID spectra were stacked together with a weight to maxi- 
mize the S/N of the final ID spectra. These ID spectra were 
calibrated to absolute flux scale by comparing them with the 
/- and //-band fluxes. The uncertainty of the absolute flux 
calibration is < 20% indicated by the difference of the scal- 
ing factor between the / and the H band. 

The cutouts from the COSMOS HST/ ACS /-band image 
( Koekemoer et al. 2007) and ID spectra of the 18 pBzKs with 
continuum detection for which the determination of the spec- 
troscopic redshifts was attempted are shown respectively in 
the left and middle panels of Figure (H 

3. SPECTROSCOPIC REDSHIFTS 

3.1. Redshift Measurement 

None of the target objects shows emission lines in the 
observed spectral range, which supports their identification 
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ID: 254025 

Zspec = 1 -8228 /zphot = 1-741 
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Zspec = 1. 4277 /zphot = 1.246 
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Zspec = 1. 4290 /zphot = 1.420 
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Figure 4. HST/ ACS stamps, MOIRCS spectra, and broad-band SEDs of the continuum detected pBzKs in this study. Left: HST/ ACS cutouts made from the 
publicly available COSMOS HST-ACS Mosaic v2.0 from Koekemoer et al. (2007). The images are 3 arcsec on a side; Middle: MOIRCS spectra. The observed 
spectrum, its 50 A median smoothed version, and its best-fit SSP template from Bruzual & Chariot (2003) library are shown respectively in gray, orange, and 
blue. The gray spectra at the bottom show the la noise. The positions of major absorption lines are indicated by the vertical lines; Right: Broad-band SEDs. The 
red circles show the observed fluxes and the gray line shows the best-fit template. The panels from the top to bottom are ordered with increasing total i^-band 
magnitude. The IDs and spectroscopic and photometric redshifts are indicated. 
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ure4. Continued. 
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ID: 269286 

Zspec = 1. 6593 /zphot = 1.665 




Figure 4. Continued. 

as PEGs. The measurement of spectroscopic redshifts was 
therefore attempted from stellar absorption features such as 
Ca II H+K, the 4000 A break and Balmer lines. The spec- 
troscopic redshifts were then derived using the Pen alized 
Pixel-Fitting method (pPXF; Cappel lari & EmsellemI 12004). 
using template stellar spectra f rom the MILES stellar library 
(ISanchez-Blazquez et al.ll2006h . We have paid special atten- 
tion to the error analysis, including systematic effects, in par- 
ticular those that could be introduced by the spurious sky 
residuals and the correlated noise. Thus, we have adopted 



the standard bootstrapping technique of resampling residuals 
to adjacent groups of spectral pixels, rather than to the indi- 
vidual pixels, to account for the correlated noise, then ran- 
domly reshuffling the best-fitting residuals in groups of 50 
pixels (300 km s~^), allowing for duplication. These residuals 
were added to the galaxy spectrum, then fitting the data with 
pPXF. To account for the sensitivity of the fits to the polyno- 
mial degree for each noise realization we adopted a random 
degree between and 4 for the additive and another random 
one for the multiplicative polynomials. In the new fits we re- 



Near-IR Spectroscopy of z > 1.4 Passive Galaxies 



9 



ID: 228121 

Zspec = 1. 8084 /zphot = 1.720 



.i (A) 




ID: 321998 

Zspec = 1-5226 / Zphot 



:t(A) 



= 1.432 




ID: 299038 

= 1.8196 /zphot = 1.767 



Arest (A) 

4500 5000 




10000 11000 



12000 13000 

Aobs (A) 



14000 15000 16000 17000 



ID: 209501 

Zspec = 1. 5409 /zphot = 1.341 




Figure 4. Continued. 

peated the entire pPXF procedure, namely redetermining the 
set of best fitting MILES stellar templates thus finding a new 
set of cr-clipped residuals. This procedure was repeated for 
300 random realizations and la errors were derived from the 
distribution of the output redshifts from the realizations. The 
results are reported in Table |3]along with their errors. 

As a cross-check, redshifts were also derived by cross- 
corre lating with simple stellar population (SSP) templates 
from lBruzual & CharlotI (l2003i hereafter BC03) with various 
ages and metallicities after broadening to a = 300 km s~^ . For 
all the objects, the resulting redshifts are consistent within la 



errors with those derived using stellar templates. 

The 18 identified objects are shown with squares in Figure 
[3l Since the identification is nearly complete (16 out of 19 
objects) at the brighter magnitude, ^ < 21, the primary factor 
of the spectroscopic identification in this kind of observations 
appears to be the total ^-band magnitude. Moreover, since we 
have selected all pBzKs brighter than ^ < 21 in the observed 
FoVs, the completeness of the spectroscopic identification is 
85% down to this depth. 

A redshift confidence class was assigned to each object ac- 
cording to the following criteria: Class 4 is assigned to the 
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ID: 253431 

Zspec = 1. 5423 /zphot = 1.558 




ID: 275414 

Zspec = 1. 8302 /zphot = 1.433 
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4500 5000 




Figure 4. Continued. 

objects which clearly show both the 4000 A break and some 
absorption lines; Class 3 is assigned to objects with a clear 
4000 A break but without unambiguously recognizable ab- 
sorption lines; Class 2 refers to redshifts derived from the 
overall shape of the continuum, i.e., the 4000 A break is not 
so prominent; Class 1 is for an insecure redshift. The assigned 
redshift classes are reported in Table |3] A redshift class > 1 
was assigned to the 18 objects for which the continuum could 
be detected, whereas for 16 objects no continuum was de- 
tected and they are listed separately in Table [T] 

3.2. Comparison between Spectroscopic and Photometric 
Redshifts 

Figure [5] shows the comparison between the spectro- 
scopic redshifts de rived above and photometric redshifts from 
lllbert et all (|2009h . This is virtually the first attempt to sys- 
tematically test the COSMOS 30-band photometric redshift 
against measured spectroscopic redshift of PEGs at z > 1 .4. A 
majority of the objects show fairly good agreement between 
photometric and spectroscopic redshifts, but several outliers 
also exist. Of course, such outliers can be ascribed to either 
of the two methods to measure redshifts, and here we briefly 
try to identify the (main) culprit. 

First notice that two outliers have quite secure spectro- 
scopic redshift, having assigned Class 4 (namely, objects 
217431 and 250093). For these objects the photometric red- 
shifts are clearly in error. In Figure |6]we show the stacked 
spectrum (see Section [5]) of the 9 objects with low confidence 
class (1 and 2): several strong features are clearly present in 
this spectrum, namely Ca II H&K, the G band and Mg/?, and 
these features are still recognizable in the stacked spectrum 
including only the five Class 1 and 2 outliers seen in Figure |5] 
We conclude that the spectroscopic redshifts are correct also 
for the majority of these two classes of objects and therefore 
in the present sample ~ 30% of the PEG photometric redshifts 
are systematically underestimated. 



ID 



254025 
217431 
307881 
233838 
277491 
313880 
250093 
263508 
269286 
240892 
205612 
251833 
228121 
321998 
299038 
209501 
253431 
275414 



Table 3 

Spectroscopic Properties 



Zspec 



Class^ 



Dn4000 



H(5f 
(A) 



1.82283 

1.42769 

1.42904 

1.8199 

1.8163 

1.4486 

1.8270 

1.52122 

1.6593 

1.54939 

1.6751 

1.42578 

1.8084 

1.52263 

1.81957 

1.5401 

1.5423 

1.8302 



±0.00055 

±0.0015 

±0.00089 

±0.0016 

±0.0038 

±0.0018 

±0.0010 

±0.00094 

±0.00058 

±0.00092 

±0.0045 

±0.00057 

±0.0015 

±0.00087 

±0.00098 

±0.0078 

±0.0037 

±0.0147 



1.40 ±0.03 
2.03 ±0.10 
2.06 ±0.08 
1.53 ±0.05 
2.74 ±0.25 
1.42 ±0.09 
1.49 ±0.10 
1.47 ±0.09 
1.38 ±0.06 
1.55±0.10 
1.52 ±0.09 
1.71±0.18 
1.38±0.12 
2.13±0.14 
2.06 ±0.12 
2.02 ±0.33 
1.61 ±0.21 
2.01 ±0.24 



2.84 ±1.30 
6.23 ±1.10 
5.59±1.15 
2.95 ±1.30 
1.30±1.51 
-2.66 ±1.78 
6.47 ±1.22 
5.40 ±1.20 
4.86±1.19 
9.03 ±0.84 
-3.92 ±1.76 
-1.80±1.69 
2.84 ±1.52 
-1.25 ±1.59 
-4.18 ±1.77 
2.07 ±1.72 
-1.47 ±1.72 
-1.59 ±1.63 



Stacked 



1.54 ±0.03 4.26 ±0.40 



^ Class 4: a spectrum with a clear detection of both absorption lines and the 4000 
A break. Class 3: a spectrum with a clear detection of the 4000 A break. Class 
2: a spectrum with a relatively high-S/N continuum detection whose overall shape 
allows redshift measurement. Class 1: a spectrum with a low-S/N continuum de- 
tection for which the derived redshift is less reliable. 

More quantitatively, the average offsets are (zspec -^phot) = 
0.12 or ((Zspec-Zphot)/(l +^spec)) = 0.04, and the standard de- 
viations of Zspec-Zphot and fepec -Zphot)/(l +^spec) are 0.14 
and 0.053, respecti vely. The normalize d median absolute 
deviation (NMAD; Hoaglin et all Il983h. (Jnmad = 1.48 x 
median(|zspec -^photl /(I +^spec)), which is used to estimate the 
accuracy of the COSMOS photometric redshift by lllbert et al.l 
(120091) . is calculated as 0.050. This is similar to ctnmad = 
0.054 f or the 1.5 < z < 3 zCOSMOS-Deep sample with 



24 (iLillv et al.ll2"QQ7h . However, they found that - 20% 
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of the 147 zCOSMOS-Deep galaxies at 1.5 < z < 3 show 
catastrophic failure defined as |zspec -^photl /(I +^spec) > 0.15, 
whereas no such large catastrophic failures are found in the 
pBzK sample presented here. Still, these r esults indicate that 
the photometric redshifts from Il bert et al.l (1200 9) would lead 
to a systematic underestimate of the volume density of high 
redshift PEGs, though the present statistics is too limited to 
precisely quantify the effect. 

Notice that all our pBzK selected objects with detected con- 
tinuum have Zspec > 1 -4, including those with Zphot < 1 -4, indi- 
cating that the passive 5z^-selection is indeed quite e ffectiv e 
to identify bona fide PEGs at z > 1.4 (see also Section [J. 2. It , 
especially when disregarding mid- and/or far-IR detected ob- 
jects. 

We also note that the clustering of the spectroscopic red- 
shifts around a few redshift spikes (including Class 1 and 2 
redshifts, see Section |4]) lends support to the correct identifi- 
cation of such redshifts. No such clustering would indeed be 
expected from random errors in the spectroscopic redshifts, 
unless spurious breaks were introduced by the flux calibra- 
tion, which does not appear to be the case. 




Figure 5. Comparison between the spectroscopic redshifts of the present 
study and the photometric redshifts from Ilbert et al. (2009). Different sym- 
bols refer to the quality Class of the derived spectroscopic redshifts as defined 
in the text and shown on the top-left comer of the figure. 



3.2.1. Improving photometric redshifts for high redshift PEGs 

To explore the possibility of improving the performance 
of photometric redshifts, specifically for high-redshift PEGs, 
we have calculated new photometric redshifts using the 18 
spectroscopic reds hifts in this work as a training sample. 
We used EAZY (Bram mer etani2008 ») to fit the 14-passband 
SEDs (UBgVrizJHKs and the 4 IRAC bands, see Section 
16.11) of the 18 PEGs, and different sets of templates, namely 
the EA ZY default templates (cf. Brammer et al. 2008, 2011; 
IWhitak er et al. 201 l|)^_aswell as passive synthetic populations 
built with Marastonrj2005L hereafter M05) and BC03 models, 
with a suita ble range of ages and metallicities. 

Following'Il bert et aH (1200 9), we estimated systematic pho- 
tometric offsets relative to the default EAZY templates by 




5000 5500 

1(A) 



Figure 6. The composite spectrum of up to 9 pBzKs whose redshifts have 
Class 1 and 2 (i.e., lower quality). The bottom panel shows the composite 
spectrum and its la noise. The locations of major absorption lines are labeled 
and indicated by dashed-lines. The top panel shows the number of objects 
used for stacking at each wavelength. 

running EAZY on the 18 PEGs having fixed the redshift to 
their spectroscopic value. This was done by measuring, for 
each band, the median offset between the observed flux of 
the galaxies and the flux of the best-fitting templates, itera- 
tively correcting all passbands until convergence. We note 
that, besides correcting for possible systematic errors in the 
photometric zero-points, aperture corrections, and/or in the 
filter transmission curves, such systematic offsets are specific 
to the adopted set of templates and reflect their possible lim- 
itations. Nonetheless, we show below that adopting the sys- 
tematic offsets determined with the EAZY templates we ob- 
tain consistent photometric redshift estimates also when using 
the passive BC03 and M05 template sets. 

In Figure |7] the resulting Zphot values are compared to the 
corresponding Zspec for the sample of our 18 PEGs, separately 
for each of the three template sets. The improvement with 
respect to the original photometric redshifts is immediately 
evident. Using EAZY default templates the median of the 
discrepancy (Zspec -^phot)/(l +^spec) is 0.006, with a scatter 
c^NMAD = 0.017. When using M05 and BC03 templates the 
median of the discrepancy is 0.004 and 0.0009, respectively 
and ctnmad — 0.02 in both cases. 

Of course, we are comparing Zspec and Zphot for the same 1 8 
galaxies used as training sample to determine the systematics 
offsets. However, within the limits due to the poor statistics, 
similarly small median discrepancies and scatters are found 
when using a random half of the sample for training and the 
other half for checking Zspec vs. Zphot- This suggests that the 
quoted accuracy of the new photometric redshifts is in fact re- 
liable, at least for sources with similar SED, magnitudes and 
redshift as the PEGs considered here. We note that the Il- 
bert et al. photometric redshifts were optimized for z ^ 1.3, 
using a set of medium bands at optical wavelengths. While 
they reached better than ~ 1% accuracy up to z 1.3 for 
bright sources, the performance for fainter sources (/ > 24) 
and par ticularly for galaxies at z ^ 1.3 is significantly de- 
graded (lllbert et al.l 12009). The larger syst ematic offset and 
scatter of the original photometric redshifts Gilbert et ani2009|) 
illustrated in Figure |5] may at least partly be ascribed to the 
lack of PEGs in the zCOSMOS-Deep spectroscopic sample 
(iLilly et al .112007 ) that was used as a training sample, as it was 
designed to contain only star-forming galaxies. 

The 18 PEGs we used for tuning photometric redshifts 
and determining their accuracy are clearly on the bright tail 
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of the general PEG population at these redshifts. Nonethe- 
less, assuming that photometric redshifts derived as above re- 
main reasonably accurate also for fainter PEGs, we measured 
the photometric redshifts for all the pBzK MIPS -undete cted 
galaxies in the COSMOS sample ("McCracken et al.'2010) and 
compared their distributio n to that obtained w ith the original 
photometric redshifts from lllbert et al.l(^009l) . 

Figure [8] shows the Zphot distribution of 90% of the full sam- 
ple, for which a reliable Zphot could be obtained as judged 
from the of the best-fit and from the Zphot probability 
distribution (dotted black and dark-gray lines, for Zphot esti- 
mated with standard EAZY and passive M05 templates, re- 
spectively). The solid lines show the distribution of the 50% 
of the sample with best constrained Zphot (i-c, with the best ^ 
values and tightest probability distribution). The dashed light- 
gray line shows th e distribution of photometric redshifts from 
lllbert et al.l (I2009h . for 85% of their full s ample (after remov- 
ing 2% of the sources unmatched in the Ilbert et al.l catalog 
or identified as X-ray sources, and a further 13% of sources 
located in masked areas). Apart from a small shift in the peak 
of the distribution (from z ~ 1.5 to z ~ 1.7), it is worth em- 
phasizing that about one th ird of the so urces are located at 
^phot < 1 .4 according to the Ilbert et al.', compared to only a 
few percent with the new determination trained on the present 
18 spectroscopic PEGs. This exercise shows that using spec- 
troscopic high-z PEGs as a training sample can substantially 
improve the estimation of photometric redshifts for this kind 
of galaxies, as illustrated by a comparison of Figure |5] and 
Figure |7l 

3.3. Note on a MIPS 24 jim source 

One of our galaxies pre-selected as a pBzK, i.e., 313880 
at z = 1.45, is detected at SpitzerMlVS 24 jim with a flux 
of 130 /iJy which corresponds to a SFR of ~ 9OM0yr~\ 
adopting a conversion based on infrared SED templates of 
IChary & Elbaz (2001) and the Chabrier initial mass function 
(IMF; Chabrier 2003). The HST/ ACS /-band (F814W filter) 
morphology in Figure |4] suggests that this object has a clumpy 
disk structure which could be the site of intense star forma- 
tion (Forster Schreiber etal.ll2011t iGenzel et al1l2011l) . Sur- 
face profile fitting (see Section [TT]) also shows this object is 
best fit by an exponential profile. Looking at the spectra in 
Figure m there seems to be an excess of flux at ~ 1.6 /im, 
indicating the presence of Ha emission as expected from the 
relatively high 24 /im flux and infrared-based SFR. However, 
in the 2D spectra taken during two nights separately as shown 
in Figure |9j there is no evidence of emission lines and the ex- 
cess is apparently caused by the residual of the OH sky line 
subtraction. Other emission lines such as [O 111] AA4959, 5007 
and H/3 are also expected in the /-band where the sky line con- 
tamination is less severe than in the if -band, but none of them 
are detected. The other explanation for non-detection of emis- 
sion lines is that the object has large amount of dust which 
attenuates the emission lines below the detection limit, (3- 
10)x 10-1'^ erg s-i cm'^ (3a). The SED fit gives Ay = 0.55!g-^4 
mag which is one of the highest among the objects in the 
present sample, but it does not seem enough to attenuate emis- 
sion lines below detection limits. Given the degeneracy of 
parameters in SED fitting, the derived Ay could possibly be 
underestimated in this particular case. Here we conclude that 
this object is probably a star-forming galaxy suffering from 
heavy dust extinction, and therefore, we excluded it from the 
stacking analysis (see Section [5]) . 



4. OVERDENSITIES OF PASSIVELY EVOLVING GALAXIES AT 
z> 1.4 

A cursory inspection of Figure |5]is sufficient to realize that 
spectroscopic redshifts are not randomly distributed, but in- 
stead tend to cluster around a few spikes. Indeed, Figure 
[TOl shows the redshift distribution of the 18 spectroscopically 
identified pBzKs in our 3 masks. There is a clear redshift 
spike of 7 pBzKs at z ^ 1.82 with la scatter of 0.006. The 
projected spatial extent of the overdensity is ~ 6 Mpc in di- 
ameter and z = 1.82 ±0.006 corresponds to 19 comoving 
Mpc along the line of sight. To our knowledge this is the first 
spectroscopically confirmed overdensity of PEGs at such high 
redshift. 

Another possible redshift- spike at z — 1.43 includes 4 
pBzKs. Although the significance does not seem high, this 
z = 1 .43 spike is also noticeable because two of its galaxies are 
the i^-band brightest pBzKs over the whole COSMOS field 
(iMancini et aLllMoh . 

A third spike appears to be present at z — 1.53, including 5 
pBzKs, and the last two galaxies have also close redshifts, 
around z — 1.67. Unfortunately, this area lies outside the 
zCOSMOS-Deep field, and therefore it is not possible to use 
its redshifts of star-forming galaxies to check for the presence 
of spikes at the same redshifts of those reported here. How- 
ever, a fair number of sBzKs was included in our MOIRCS 
masks (cf. Table [2]) and a comparison of the redshift distribu- 
tions of pBzKs and sBzKs will be done in a future paper. 

Figure [TT] shows the location on the sky of the galaxies in 
these three redshift spikes at z = 1 .43, 1 .53, and 1 .82, together 
with that of other photo-z selected galaxies with Az = ±0.02 
and ^ < 23 from the mean redshift of each spike. There may 
well be large spikes at these redshifts extending spatially well 
beyond the area explored by the three MOIRCS FoVs, but our 
data do not allow us to claim the presence o f bound clu sters. 

The XMM-Newton data (Hasinger et a l.l 120071; 
iFinoguenov et al.l [20071) indicate that there is no extended 
X-ray emission associated to these redshift spikes down to 
(6-8) X 10-15 erg s'^ cm'^ (2a) in 0.5-2 keV band, which 
corresponds to the 2a limit of Lx — (1-2) x 10"^^ erg s"^ 
and M200 ^ 1.5 X lO^^^M© at the redshifts of the spikes 
(A. Finoguenov, private communication). 

5. THE COMPOSITE SPECTRUM OF pBzKs 

A composite spectrum has been constructed by stacking 
17 out of the 18 spectroscopically identified pBzKs. We ex- 
cluded the object 313880 because of the MIPS 24 /im detec- 
tion (see Section [33]) , which could cause a contamination of a 
star-forming galaxy into a passively evolving population. We 
did not exclude objects with Class 1 and 2 since the spec- 
troscopic identifications appear to be correct as investigated 
in Section 13.21 and low-S/N does not affect the result after 
weighting as described below. 

After being de-redshifted to the rest-frame, each individual 
spectrum was first normalized by the mean flux at rest- frame 
4500 < A(A) < 5200. Then the spectra were linearly inter- 
polated into a 1 A linearly spaced wavelength grid. The as- 
sociated noise spectra were normalized by the same factor as 
that for the object spectrum and interpolated to the rest- frame 
in quadrature. The spectra were then coadded with weights 
proportional to the inverse square of the noise spectra at each 
rest- frame wavelength. 

To estimate a realistic error and any biases related to the 
stacking, we applied the jackknife method to the spectra used 
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Figure 7. Comparison between the spectroscopic redshifts from the present study and the corresponding photometric redshifts from different determinations: 
Left: Zphot fr om this work using standard EAZY templates; Middle and Right: Zphot from this work using EAZY with passive BC03 and M05 models, respectively 
(see Section [3.2. IV All EAZY photometric redshifts were obtained using the 18 PEG spectroscopic redshifts in the present study as a training sample. The 
solid line traces the bisector whereas the filled areas show deviations of 0.025, 0.05, and 0.1 x (1 Readers can refer to Figure |5]for a comparison with the 
photometric redshift by Ilbert et all (120091) . 
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Figure 8. The distribution of photometric redshifts for pBzK-selected PEGs. 
The solid/dotted lines show the redshift distribution of 50%/90% of the full 
PEG sample over the COSMOS field using the photometric redshifts mea- 
sured in this work using different sets of templates (standard EAZY and pas- 
sive M05 templates in black and dark gray, respectively). The dashed light- 
gray line shows the redshift distribution of 85% of the full PEG sample using 
flbert et al. (2009) photometric redshifts. All EAZY photometric redshifts 
were obtained using the 18 PEG spectroscopic redshifts in this paper as a 
training sample. 




Figure 10. The redshift distribution of spectroscopically identified pBzKs 
in the present sample. 
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Figure 9. The two dimensional spectra of the MIPS 24 fim detected source, 
313880, taken on February 6, 2010 (Top) and on April 1, 2010 {Bottom). The 
spectra are 1.9 arcsec wide in the spatial direction. 

for stacking. We made 17 composite spectra in the same way, 
but removing one object at a time from the stacking. Then 
the standard deviation of the flux at each wavelength pixel is 
estimated as 

-Jack ~E (/-/(,))', (1) 



where A/^ is a number of objects in the sample (i.e., N = 17 
here), / is a flux of composite spectrum of N spectra, and /(/) 
is the flux of the composite spectrum made oi(N-l) spectra 
by removing /-th spectrum. We use ajack as the la error of 
the composite spectrum. The composite spectrum was also 
corrected for a sampling bias derived from the jackknife esti- 
mator as 

/ = /-(yv-i) ((/(,))-/), (2) 

where (/(/)) is the average of /(/). The typical correction fac- 
tor due to the bias estimate is < 2%. We adopt the bias cor- 
rected spectrum f as the final composite spectrum, which is 
inevitably dominated by the objects with highest S/N ratios as 
well as by observed-frame wavelength less contaminated by 
the sky residuals. 
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Figure 11. The spatial distribution of the pBzKs in the three major redshift 
spikes. The symbols are redshift-coded as indicated in the insert. Small 
blue pluses, green crosses, and orange circles show objects with K <23 and 
photometric redshift within Az = ±0.02 from the mean redshifts of the spikes 
at z = 1.43, 1.54, and 1.82, respectively. Solid boxes and circles represent the 
position of MOIRCS detectors and masks, respectively. 



The composite spectrum obtained in this way (equivalent 
to a MOIRCS integration time of about 140 hours) is shown 
in Figure [121 At least 1 1 spectral lines/features are clearly 
visible in this spectrum, and identified in Figure [121 specifi- 
cally, CN, Ca II H&K, m, the G band, H7, Fe4383, Fe5270, 
Fe5335, H/3, and Mgb. Although the number of stacked 
objects are small at longer wavelength where contamination 
from the residual OH-airglow subtraction is severe, there 
could be an indication of excess flux at the location of Ha. 
We have measured the equivalent width (EW) of the poten- 
tial Ha emission line from th e sta cked spectrum divided by 
the best-fit model (see Section [64l) accounting the underlying 
absorption, and we found EW(Ha) = 3.5-5 A depending on 
the width of the integration window to derive the EW (1-2.5 
times of FWHM of the instrument). In case of EW = 5 A, it 
can be translated into SFR = 5.6 Mq yr~^ and 1 .9Mq yr~^ for a 
object at z = 1 .69 with the brightest and median H-band mag- 
nitudes, respectively. In addition to Ha, we have also carried 
out a similar measurement on the position of [O II] A3727, and 
found corresponding SFR of 4M0 yr~^ and 1 .6Mq yr~^ for the 
brightest and median cases with an error of about 20-40%. 
Looking at the individual spectra extending to Ha, 307881 is 
the brightest and could mainly contribute to this feature in the 
composite spectrum. However, the feature around Ha in the 
spectra of 307881 appears to be rather broad, mostly contami- 
nated by adjacent OH residuals. Indeed, when we remove this 
object from the stacking, the feature disappeared, while the 
estimated SFR at the position of [O II] still remains the same. 
Therefore, the residual star formation, if any, could be at most 
~ 5M0yr~^ and most likely < l-2M0yr~\ which is > 100- 
200 times smaller than those of galaxies with = 10^^ M^:^ 
on the star formation main sequence at z = 2 (e.g.,lDaddi et al.l 



120071: iPanneUa et al.ll2009l:lRodighiero et al.ll2011h . The com- 
posite spectrum is analyzed is Section [631 



6. STELLAR POPULATION PROPERTIES OF THE PROGRAM 
GALAXIES 

6. 1 . Broad-band SED Fitting 

Having fixed the redshifts to their spectroscopic values, we 
have derived the physical properties of the galaxies from their 
broad-band SEDs. Besides the Subaru/Suprime-Cam B and 
z ' band and CFHTAVIRCam ^-band photometry to apply the 
BzK selection technique, we used the multi-band photometry 
including the CFHT/MegaCam u band, the Subaru/Suprime- 
Cam g\ y, r\ and /' bands, and the CFHTAVIRCam / and 
H bands as well as the four Spitzer/IRAC bands. All mag- 
nitudes from uio K are converted to total magnitudes by ap- 
plying aperture corrections both for point sources and for ex- 
tended objects (McCracken et al. 2010; Mancini et al. 2011). 
The IRAC band magnitud es are c orrected to the total magni- 
tude following Il bert et all (l2QQ9l) . 

The SED fitting was carried out b y using FAST (Fitt ing and 
Assessment of Synthetic Templates : [Kriek et al.ll2009l) . FAST 
includes internal dust reddening and provides various physical 
quantities of the galaxies such as stellar mass, age, and metal- 
licity. For the templates, we used the composite stellar popu- 
lation models generated from the SSPs of Chariot & Bruzual 
(2007, in preparation; hereafter CB07) with the Chabrier IMF. 
It has been argued that the IMF may not be universal among 
local early-type galaxies, but could depend on the galaxy stel- 
lar mass (Cappellari et al. 2012). If so, the masses derived 
here would be underestimated by up to a factor 2-3 for the 
most massive objects. We assume a universal IMF in this pa- 
per, as it was done in all other similar studies to which we 
compare our results. 

The adopted star formation histories (SFHs) assume an 
exponentially declining SFR with various ^-folding times 
(r) ranging from log(r/yr) = 7.0 to 10 with an interval of 
Alogr = 0.1 and ages ranging between log(age/yr) = 7.0 to 
10 with an interval of A log age = 0.05. The metallicity Z is 
fixed for each model (no attempt is made at mimicking chem- 
ical evolution) and is chosen from the set Z = 0.004, 0.008, 
0.02, and 0. 05. Dust extinction was also applied following 
the recipe bv lCalzetti et al.1 (l2000h with Ay = 0-3 mag with an 
intervals of 0.05 mag. Ages are required to be less than the 
age of the Universe at the observed redshift. The time-like 
free parameters in the fit are therefore the age and r. 

Then the best values and 68% confidence intervals of the 
stellar mass, timescale of the star-formation (r), age from the 
onset of star formation, reddening, and metallicity for each 
galaxy are derived based on the likelihood distribution of each 
template. The SFR- weighted ages are also computed by using 
the best-fit r and age. Derived physical parameters are listed 
in Table[4land the observed galaxy SEDs and the best-fit tem- 
plates are shown in the right panels of Figured 

All but a few of the pBzKs have SEDs consistent with 
Ay < 0.4 or E(B-V) < 0.1, i.e., as expected the derived dust 
extinction is very small for most of the objects. When con- 
sidering that reddening and age are partly degenerate one can 
conclude that none of our objects suffers major dust obscu- 
ration, as expected for passive galaxies. Stellar population 
age of the pBzKs are typically > 2 Gyr. All best-fit SEDs 
have a short ^-folding time r < 500 Myr which is typically 
much shorter than the derived age, i.e., the present star for- 
mation is already well suppressed. Stellar masses range from 
^ 4 X 10^^ Mq to ~ 4 X 10^^ M0, whereas best fit metallici- 
ties are near- solar or slightly sub- solar with large uncertainties 
that cover well above the solar metallicity. 
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Figure 12. The same as Figure|6] but now co-adding the spectra of all 17 MIPS -undetected pBzKs in the present sample. 



It has been recently argued that an increasing SFR is more 
appropriate than a decHning one in the case of high-redshift 
star- f o rming galaxies (e.g., iRenzinil I2QQ9I: [Maraston et al.l 
120101: iPapovich et al.l |2011|) . In particular, an exponen- 
tially increasing SFR has been suggested, with SFR cx 
r~^exp(+^/r) where r 500 Myr. In the case of pBzKs 
such SFR needs to be truncated at some point, in order mimic 
the quenching of their star formation. Most galaxies in the 
present sample show a strong 4000 A break a nd so me of them 
strong B aimer absorption lines (see Section 16. 3K indicating 
that the age of the dominant stellar populations must be > 1 
Gyr, hence star formation has to be quenched > 1 Gyr before 
the observed epoch. As shown in Table HI assuming instead a 
declining SFR the best fit star formation timescale turns out to 
be typically about 10 times shorter than the age of the galax- 
ies which is well in excess of 1 Gyr. Therefore, for these 
quenched galaxies both exponentially declining and increas- 
ing star formation scenarios would result in almost identical 
stellar population properties being selected by the best fit pro- 
cedure. 

6.2. Rest-frame UV J colors 

Passively evolving galaxies are known to lie in a distinct 
region on rest- frame ( U -V) vs. (V- 7) two -color diagram 
(e.g., Wu vts etal.1 120071: IWhitaker et ani2011h . This can be 
used to check the passive nature of pBzKs in the present sam- 
ple. Figure[T3]shows the UVJ diagram where we have adopted 
the identical filter combination as that used in'Whitaker et al.' 
fcOll), namely the Johnson UBV system (Maiz Apellaniz 
[2006) and the 2MASS / band. The rest-frame UVJ colors 
were computed by convolving the best-fit templates with the 
filter response curves without applying atmospheric and in- 
strumental throughputs. The selection box for passive galax- 
ies is defined as {U-V) > 0.8 x (y-/) + 0.7, (U-V) > 1.3, 
and (V-J) < 1.5, and young and old galaxies are further sep- 
arated at (V-J) = 0.9 following .WhitakeretaL (-2012) . 



As shown in Figure[T3l most of objects have rest- frame UVJ 
colors consistent with passive evolution. The colors of some 
outliers are still within '^0.2 mag from the selection box, and 
have an offset consistent with possible systematics in the SED 
fitting and photometric errors. In fact, galaxies with stellar 
ages of > 1 Gyr can still be found around (U-V) 1.7 and 
(V-J) 1.5 (Whita ker et al.ll2012) . The MIPS 24 /imsource, 
313880, is well within the selection box for passive galaxies, 
which suggests that the object may be a contaminant to the 
UVJ selection, or is affected by systematics in the SED fitting. 
From this UVJ color test we conclude that in all 1 8 objects star 
formation has been already quenched, or is very low. 

6.3. Rest-frame Optical Absorption Line Indices and the 
Dn4000 Index 

To further characterize the stellar population content of 
the pBzKs in this sample some of the Lick indices and the 
strength of the 4000 A break were measured on our rest- 
frame optical spectra. To this end we used the Lick_EW 
program as a part of the EZ_Aqes IDL code packag43 
(iGraves & SchiavonI 120081: [Schiavon 2007), following the 
definition o f the Lic k indic es by IWorthey et al.l (Il994|) and 
f^rthe y & Ottavianil (|1997|) . For this purpose one needs to 
specify a stellar velocity dispersion and the instrument disper- 
sion: we adopted 27 A in the observed frame for the FWHM 
of the instrument configuration, and a stellar velocity dis- 
persion of cr* = 270 km s~^ for 254025 (see Section 17.41) . 
For all other objects we used the stellar mass-stellar veloc- 
ity dispersion relation of local galaxies at 0.05 < z < 0.07, 
determined by using stellar velocity dispersions from the 
NYU Value-Added Galaxy Catalog (Blanton et al. 2005) and 
stella r masses from the MPA- JHU catalog ( Kauffm ann et all 
120031: ISalim et al.ll2007h . The stellar masses of SDSS galax- 
ies are converted to a Chabrier IMF from a Kroupa IMF 

^ Jhttp: //astro .berkeley . edu/~ { } graves/ez_ages .html] 
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Table 4 

Physical Parameters from Broad-band SED Fitting 
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254025 
217431 
307881 
233838 
277491 
313880 
250093 
263508 
269286 
240892 
205612 
251833 
228121 
321998 
299038 
209501 
253431 
275414 
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1 1 05+0-08 
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^-0.24 




008+0 018 
0.008f 
0.008+ 
0.008+ 
0.020+ 
0.008+ 
0.020+ 
0.020+ 
0.004^ 
0.020+ 
0.020+ 
0.020+ 
0.020+ 
0.020+ 
0.008! 
0.020+ 
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^ Age from the onset of star-formation. 

^ Age weighted by SFR derived as fQ^{tQ-t)ip{t)dt / ip{t)dt, where ip{t) is the SFR and to 
and r are the best-fit age and r, respectively. 
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Figure 13. The rest-frame ( U - V)-{y - J) diagram for the 1 8 pBzKs in the 
sample (circles). The colored areas represent the selection box for young 
(blue) and old (red) passively evolving galaxies, respectively, following the 
definition bv iWhitaker et al. (i2012i) . 



iKroupal (1200 ll) , albeit the correction factor for the veloc- 
ity dispersions is small (~ 10% and 3.5% increases for the 
stellar mass and stellar velocity dispersion, respectively). It 
is quite possible that high redshift PEGs are more compact 
than local ones (cf. Section 17.11 below), hence may have 
higher velocity dispersion compared to the local M^-a re- 
lation. However, we use here the local relation because 
with the exception of the value reported by ivan Dokkum et al.l 
(I2OO9.) . all other objects at z > 1.4 with measured a still fol- 
low such relation (C appellari et al. 2009; Cenarro & Trujillol 
120091: iQnodera et al.ii20m:.van de Sande et al. 2011). 

Given the low-S/N of our spectra (typically S/N < 5 in the 



60 km s~^ interval) and the pixel to pixel fluctuations in S/N 
due to sky emission lines, we modified the Lick_EW code 
to use median flux values to compute the pseudo-continua for 
the indices. Moreover, absorption lines are not clearly visible 
in most of the spectra, and therefore we restricted our analysis 
only to objects in which lines are clearly detected and allow 
a fairly accurate measurement of the relative index. Measure- 
ments have been carried out for all objects. In particular, here 
we focus on the H(5f index, which is less affected by OH lines 
and by emission line filling, if any. Ideally, H/3 would be a 
more useful index because it is relatively isolated from other 
metal absorption lines, hence independent of metallicity of 
the stellar population. However, for many of our objects H/3 
is strongly contaminated by sky emission lines. 

Besides the H(5f index, we measured the strength of 4000 
A break as quantified by the Dn4000 inde x, for which we 
follow the definition of iBalogh et al.l (119991) which uses rel- 
atively narrow wavelength windows (A A = 100 A) for red and 
blue continua. The Dn4000 index is more robustly derived 
because the wavelength intervals are wider than those used 
for absorption line indices. Errors on Dn4000 were estimated 
through Monte Carlo simulations by artificially adding ran- 
dom noise assuming the normal distribution based on the error 
spectra at each wavelength pixel. Then 6 8 -percentile intervals 
are derived from 10,000 realizations per object. 

Figure [14] shows the zoom-in of the spectra around H(5 for 
the objects with high-S/N (S/N > 5) for which the absorption 
line is relatively well detected. On the other hand. Figure [15] 
shows the spectra around the 4000 A break for all objects with 
spectroscopic redshift. For almost all the objects, Dn4000 ap- 
pears to be measured reasonably well, while for the objects 
not shown in Figure[T4]the region around H^ is dominated by 
the noise. 

Figure [16] shows the Dn4000 and Yid^ as a function of the 
stellar mass, along with the corresponding relation for lo- 
cal galax ies at z 0.06 extracted from the MPA-JHU SDSS 
database (iKauffmann et al.|[2003h . corrected for emission line 
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filling. Local galaxies fall in two sequences, one for passive 
galaxies with higher Dn4000 and lower HJp at a given stel- 
lar mass, and another for star-forming galaxies with lower 
Dn4000and higher H(5f. 

The objects with S/N > 5 are also the most massive pBzKs 
in the sample. As illustrated in Figure [17] which shows HS^ as 
a function of Dn4000, two among them (217431 and 307881) 
show simultaneously high Dn4000 (~ 2, comparable to those 
of the local early-type galaxies) and high H(5f 5 A), much 
stronger than typical of local PEGs of similar stellar mass. 
A combination of high Dn4000 (> 1.5) and high 115^ (> 3 
A) appears to be typica l of local and mode rate redshift post- 
starburst galaxies (e.g., Vergani et al.'"2010j), although in this 
particular case the combination (Dn4000, H(5f) — (2, 5) ap- 
pears to be rather extreme (cf. e.g., Le Borgne et al. 2006). 
This can be also seen in Figure [17] in which the location of 
these two objects cannot be explained by any of the overplot- 
ted star formation histories, namely instantaneous burst, expo- 
nentially declining SFR, constant and exponentially increas- 
ing SFR with truncation. We cannot exclude that Dn4000 may 
have been somewhat overestimated, given the relatively nar- 
row wavelength range available shortwards of the break (see 
the corresponding spectra in Figure |4), although the sensitiv- 
ity curves used for the flux calibration are still smooth and 
do not suffer from a sharp drop at the edge. By taking these 
indices at face value, we suggest that these two galaxies may 
have been quenched a relatively long time ago compared to 
the star formation timescale, which makes Dn4000 larger, and 
experienced a recent episode of star formation which has en- 
hanced the absorption, although it appears to be difficult 
to enhance H(5f without changing significantly Dn4000. 

The spectra of the two other objects with high S/N and the 
composite spectrum show a Dn4000 and HJf combination 
more akin that of local star-forming galaxies. However, there 
are virtually no local star forming galaxies with comparable 
stellar mass and Dn4000 as at this stellar mass almost all local 
galaxies are already passive since many Gyr. This is also the 
case for HSj^ as there is almost no local counterpart with sim- 
ilar stellar mass and H(^f. The location of these two galaxies 
in Figure [17] i.e., Dn4000 1.5 and H(5f A, are typical 
of the early stages of a post-starburs t, i.e., of a very recently 
quenched galaxy (iBalogh et al.ll 19991) . 

6.4. Stellar Populations and Star Formation History from the 
Composite Spectrum 

The high S/N of the composite spectrum allows us to study 
the "average" stellar population content of our PEGs at z > 1 .4 
by using both the detailed shape of the continuum and the 
absorption lines. For this purpo se, we used STARLIGHT0 
(ICid Femandes et al.ll2QQll2009l) . STARLIGHT is a program 
to fit an observed spectrum with a combination of template 
spectra which usually consist of population synthesis mod- 
els. These template spectra are selected from the CB07 SSP 
library spanning 29 ages logarithmically spaced between 10 
Myrs and 5 Gyrs (the age of the Universe at z = 1.4 is ~ 4.5 
Gyrs) and 4 metalHcities from Z = 0.004 to Z = 0.05. No in- 
trinsic reddening was assumed as there are no emission lines 
in the composite spectrum indicative of a star formation activ- 
ity. We leave the total velocity dispersion and velocity offset 
as free parameters. 

The results of the fitting by STARLIGHT are shown in Fig- 

^ |http : / / www . starlight . uf sc . bT] 



ure [18] The composite spectrum is well reproduced by the 
model with the rms of the residual being < 8% of the observed 
flux and with a reduced-x^ = 1.02. The composition of stellar 
populations giving the best fitting spectra is shown in the left- 
most panels of Figure[T9] The light at 4020 A is dominated by 
400-800 Myr stellar populations with Z = 2.5Z0, whereas the 
mass is dominated by 3 Gyr old stellar populations, again 
with Z = 2.5Z0. The luminosity and mass weighted ages are 
(log^/yr)^ = 9.01 and (log^/yr)^ = 9.59, respectively. 

To check the dependence on the choice of template spectra, 
the composite spectrum was fitted by restricting the metal- 
licity within the range Z = 0.008-0.02 and/or dropping the 
templates with age less than 100 Myrs. The quality of the fit 
is essentially indistinguishable from the previous one, with 
a reduced-x^ = 1.02-1.14 and the resulting SFHs are also 
shown in Figure [19] There is no indication for the pres- 
ence of young stellar population with age of < 400 Myr in 
all fitting results, i.e., there is no detectable B-type star con- 
tribution, consistent with the strong B aimer absorption lines 
(mainly produced by A-type stars) seen in the composite spec- 
trum. The lack of B-type stars is also consistent with no emis- 
sion line being detected in the residual spectrum in Figure [18] 
When the metallicity range is limited to sub- solar (two right- 
most panels, using only templates with Z = 0.004 and 0.008), 
the composite spectrum is reproduced by a ~ 1 Gyr old stel- 
lar population formed in an almost instantaneous burst. On 
the other hand, if solar or super-solar metallicities are allowed 
in the fitting, the mass results to be dominated by > 3 Gyr old 
stellar populations with a small amount of younger, ~ 400- 
800 Myr old stellar populations, very similar to SFR- weighted 
ages derived from SED fitting. 

An independent check has been carried out by 
doing the spectral fitting using the pPXF routine 
(ICappellari & EmsellemI 120041) . T he fit used as tem- 
plates a set of 210 MIL ES SSP JVazdekis et alJIgOlOj: 
Falcon ^arroso et aT] 1201 Ih with a Kroupa IMF iKroupal 
(2001) arranged in a regular grid of 35 different ages, loga- 
rithmically spaced between 0.1 and 5 Gyr, and 6 metallicities, 
namely [M/H] =-1.71, -1.31, -0.71, -0.40, 0.00, and 0.20. 
To reduce the noise in the recovery of the SFH and metallicity 
distribution we empl oyed linear regularization of the weights 
during the fitting ( Press et aH [1992), as implemented in 
the current version of pPXF. The regularization level was 
adjusted to increase the x^ from the best un-regularized fit 
by Ax^ = Y^2 X A^pix, with A/pix the number of spectral pixels. 
In this way the solution represents the smoothest one that 
is still consistent with the observed spectrum. The obtained 
result is shown in Figure [20] as a mass fraction within each 
age and metallicity interval. Since the best-fit spectrum is 
indistinguishable from that from STARLIGHT, we do not 
show it here. The inferred distribution is dominated by 
a ~ 2.5 Gyr stellar populations with solar to super- solar 
metallicities, with at most a small (< 1%) contribution from 
young (100-200 Myr) populations. The mass- weighted age 
is (log^/yr)M = 9.40, which is very close to that found using 
STARLIGHT. 

These results show once more the effects of the well known 
age-metallicity degeneracy, in particular, from the fit with 
pPXF which produces nearly constant contours along a line 
with decreasing age and increasing metallicity. Combining 
the two completely independent experiments, derived by us- 
ing different codes and spectral libraries, any significant con- 
tribution from young populations can be excluded and we can 
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Figure 14. Zoom-in of the spectra around 116 for the four objects with the highest S/N ratio. The soHd Hne in each panel shows the object spectrum and the 
corresponding la noise. The left and right gray-shaded regions indicate respectively the blue and red continuum windows used used in the definition of the HSf 
index, whereas the central ones show the window including the ¥16 absorption line. The horizontal thick lines indicate the level of blue and red continua for each 
object. The noise spectrum is shown at the bottom of each panel. 
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Figure 15. Zoom-in of the spectra around the 4000A break. The solid line in each panel shows the spectrum and their la noise is shown by the light-yellow 
shaded regions. The gray shaded regions indicate the blue and red continuum windows used to measure the Dn4000 index, and the horizontal thick lines indicate 
the adopted continuum level on each side of the break. The brown shaded region at the bottom of each panel represents the 1 a noise divided by VlOO (the size 
of each continuum window), which shows the level of the effective noise associated to the flux in these continuum windows. 

conclude that the composite spectrum is well fitted by old (> 1 
Gyr), metal-rich stellar populations. The very small contribu- 
tions by < 400 Myr old populations found by STARLIGHT 
and by ~ 200 Myr old populations found by pPXF are not 
significant, and should rather be interpreted as upper limits. 



7. STRUCTURAL AND KINEMATIC AL PARAMETERS 

A puzzling property of high-redshift PEGs has emerged 
in recent years and is currently much debated in the lit- 
erature. Indeed, several of them appear to have effective 



radii ^ 2-5 smaller compared to local ellipticals of the 
same st ellar m ass (e.R.,'Daddi et al. "20051; iTruiillo et alJ[2007l; 
Cimat ti et al.l [20 08: van Dokkum et al. 2008; Sara cco et aO 
i2009i: ICassata et al.. 2011), which implies densities within 
their respective effective radius that are ~ 10-100 times 
higher than that to their counterparts at z = 0. No general 
consensus has yet emerged about the physical mechanisms 
that would cause such observed size evolution, although ac- 
cretion of an extended envelope via minor mergers is one 
widely-proposed scenario. On the other hand, z > 1.4 PEGs 




Figure 16. Left: the Dn4000 index as a function of stellar mass for the 18 pBzKs with spectroscopic redshifts. Right: the H(5f index as a function of stellar 
mass. The number of objects in each bin of Dn4000 and H(5f are shown as histograms on the right panel of each plot. Individual objects are shown with red 
circles, yellow triangles, and white squares corresponding to S/N > 5, 3 < S/N < 5, and S/N < 3 in 60 km s~^ bins, respectively. The pink pentagon refers to the 
composite spectrum where the stellar mass is a S/N-weighted average among those used for the stacking. The underlying gray shaded area shows the distribution 
of the local SDSS galaxies including both passive and star-forming objects, after correction for emission line filling in the case of emission line galaxies. 
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ture 17. The YiSp vs. Dn4000 plot; symbols correspond to those in Figure 
The contours show the distribution of the local SDSS passive galaxies 
whereas the gray shaded area refers to both passive and star- forming local 
SDSS galaxies. The solid lines show the evolution of stellar populations with 
various star formation histories (SFH): an instantaneous burst (i.e., a SSP, 
blue), exponentially declining SFHs with ^-folding time of 1 Gyr (purple), 
constant star formation rate over 1 Gyr followed by instantaneous quenching 
to SFR = (green), and exponentially increasing star formation rate with 
T = 0.5 Gyrs with a truncation at 1 Gyr (orange), as also indicated in the 
insert. The CB07 evolutionary synthesis models have been used. 

wit h size comparable to that of local ellipti cals do also ex- 
ist (|Mancini et al.ll201Ql: ISaracco et al.ll201Ql) . indicating that 
galaxies with a diversity of structural properties co-exists at 
these high redshifts. Possible biases that could lead to size 
underestimates have also been discussed (e.g ., iDaddi et aH 
|2005; Hopkins et al. 2009; Mancini et al. 201Qh. but it is un- 
derstood that observational biases could account for only part 
of the effect. An independent way to check this issues is 
by measuring the stellar velocity dispersion (a*) from ab- 
sorption lines. If they are truly compact and highly dense, 
their cr* should be higher than that of local ellipticals. So 




Figure 18. The STARLIGHT fit to the composite spectrum. The gray line 
is the composite spectrum being thicker in the wavelength range used for the 
fit. The orange line shows the STARLIGHT best fit model having allowed 
age and metalHcity to span the 0.01 < Age(Gyr) < 5 and 0.004 < Z < 0.05 
ranges. The gray line around zero gives the difference between the model and 
composite spectrum. 

far cr* has been measur ed for onl y 5 or 6 ind ividual PEGs 
at z > 1.4 (CapDellari efaT'IOO^'; Cenarro & Truiillo 2Q0l 
2009; 



van Dok kum et al. 2009; Newman et a l. 2010; Onodera e t aL. 
20101 van de Sande et al.l l20TTl) . and with one exception 



(Ivan Do kkum et al.. .2009) they appear structurally and dy- 
namically similar to local ellipticals. 

The structural parameters of the galaxies in the present sam- 
ple are then presented and discussed in this section. We at- 
tempted the extraction of the velocity dispersion with pPXF 
for all galaxies. However, due to (i) the low S/N, combined 
with (ii) the presence of systematic sky residuals, and (iii) the 
low spectral resolution, a reliable results could be obtained 
only for the galaxy with the highest S/N. 

7.1. Surface Brightness Fits and Measurement of Galaxy 
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Figure 19. The age and metallicity distributions of the stellar populations resulting 
refer to the fractional contributions to the light at 4020 A (top) and mass (bottom). 
indicated, whereas metallicities are color-coded on the bottom-left panel. 

Mass Fraction 



0.0 



^ -1.0 
-1.5 



0.05 
0.04 
0.03 

0.01 
0.00 



0.1 



1.0 

Age (Gyr) 



Figure 20. 

fitting to th 

The colors indicate the mass fraction in each age and metallicity bin 



The age and metallicity distributi on obtained from the spectra l 
fitting to the composite spectrum with pPXF (Cappe llari & Emsellemll2004r) . 



We have measured the structural properties of these galax- 
ies by fitting the 2D Sers ic profile with GALFIT version 3.0 
(iPeng et al.ll2QQl IMOah on the HST/ACS /-band (F814W) 
mosaic version 2.0 over the COSMOS fi eld which has a 
reduced pi xel sca le of 0.03 arcsec pix~^ (iKoekemoer et al.l 
[2007; Ma ssey et al. 2010), obtained f rom the original i mages 
with the MultiDrizzle routine (K oekemoer et al.ll2002|) . The 
structural properties of three of these 18 galaxies (namely 
25 4025, 307881 and 217431) have already been measured 
bv lMancini etaH dlOlO.) on the previous COSMOS HST/ACS 



I from the STARLIGHT best fit to the composite spectrum. The distributions 
On each top panel the age and metallicity ranges allowed during the fit are 

release which used the original ACS pixel scale of 0.05 arc- 
sec pix~^ . This difference has slightly affected the GALFIT 
results (see the next section). Here we have paid special atten- 
tion to the treatment of neighboring objects, such as in object 
254025 which has one bright and two faint neighbors. There- 
fore, in this and similar cases we used the segmentation map 
generated by SExtractor (Berti n & Amo uts 1996 ) to mask out 
the faint neighbors and fit the brightest ones together with the 
main object. 

For each galaxy, we constructed the PSF by combining un- 
saturated stars in the field as close as possible to each target, 
though using the single nearest star does not change the re- 
sult at all. Indeed, the Sersic index n and effective radius re 
vary on average within 2% of the values reported in Table |5] 
Some of the objects are not bright enough in the /-band (e.g., 
/ab > 24) to allow a robust estimate of the structural param- 
eters, though we have carried out the measurement for the 
whole sample. 

The level of the sky background is the most crucial pa- 
rameter which affects the estimates of rg and the Sersic in- 
dex whereas the estimate of the /-band magnitude, position 
angle, and axial ratio are less affected by the adopted sky 
level. To account for the effect of the sky background, we 
ran GALFIT with different assumptions for the sky value: (i) 
sky as free parameter in the fit, (ii) sky fixed to the so-called 
"pedestal" GALFIT estimate, and (iii) sky manually measured 
from empty regions near the main object. By combining the 
results from such GALFIT runs, we derived rg, n, /-band mag- 
nitude, and axial ratio as the mid point between their maxi- 
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Table 5 

Structural Properties 











Sersic n 


^814 




(arcsec) 


(KpC) 




(mag) 




(2) 








P) 


254025 


0.375 ± 


0.072 


3.16±0.61 


3.36 ±0.42 


22.83±0.13 


217431 


0.851 ± 


0.231 


7.19±1.95 


3.83 ±0.56 


22.38±0.18 


307881 


0.318± 


0.014 


2.68±0.12 


2.29 ±0.10 


22.92 ±0.04 


233838 


0.266 ± 


0.037 


2.25 ±0.31 


3.07 ±0.34 


23.74 ±0.09 


277491 


0.292 ± 


0.013 


2.46 ±0.11 


1.00^ 


24.31 ±0.05 


313880 


0.460 ± 


0.063 


3.89 ±0.53 


1.00^ 


23.53±0.19 


250093 


0.356 ± 


0.018 


3.00±0.15 


1.00^ 


24.33 ±0.05 


263508 


0.101 ± 


0.004 


0.86 ±0.03 


3.23 ±0.21 


23.76 ±0.03 


269286 


0.122± 


0.014 


1.03±0.12 


4.96 ±0.72 


23.82 ±0.06 


240892 


0.153 ± 


0.016 


1.29±0.13 


3.01 ±0.34 


24.24 ±0.07 


205612 


0.245 ± 


0.033 


2.08 ±0.28 


2.36 ±0.26 


24.17±0.11 


251833 


0.190± 


0.011 


1.61 ±0.09 


2.13±0.13 


24.01 ±0.04 


228121 


0.329 ± 


0.102 


2.78 ±0.86 


4.07 ±0.89 


24.40 ±0.20 


321998 


0.216 ± 


0.040 


1.83 ±0.34 


4.40 ±0.65 


24.15±0.10 


299038 


0.114± 


0.002 


0.96 ±0.02 


1.89 ±0.05 


24.94 ±0.02 


209501 


0.077 ± 


0.000 


0.65 ±0.00 


1.63 ±0.01 


24.97 ±0.00 


253431 


0.317± 


0.022 


2.69±0.18 


4.23±0.19 


24.46 ±0.04 


275414 


0.269 ± 


0.015 


2.27±0.13 


1.00^ 


25.07 ±0.05 



Note. — (1) ID; (2) effective radius in arcseconds; (3) effective radius in kpc; 
(4) Sersic index; (5) HSTIACS /-band magnitude recovered by GALFIT. The effective 
radii are circularized and the disk scale lengths for the objects best fitted with the 
exponential profile are converted into effective radii. 



^ The fitting has been carried out by fixing « = 1 for these objects. 

mum and minimum values and quoted half such range as the 
corresponding error. This method provides a more reliable e s- 
timate of the actual uncertainties of these measurements com- 
pared to the small formal GALFIT errors from the test. 

Four of the 18 galaxies (namely, 250093, 275414, 277491, 
and 313880) are better fit with a pure exponential profile 
rather than a Sersic profile with n> I. For these objects the 
normal GALFIT fit did not converge when leaving ^ as a free 
parameter, whereas an acceptable fit can be achieved when n 
is fixed to unity. The resulting structural parameters are re- 
ported in Table [5j where effective radii are circularized (i.e., 
Te = \lhja where is the effective semi-major axis as given 
by GALFIT and is the axis ratio). The disk scale lengths 
for the objects best-fit with the exponential profile are trans- 
lated into effective radii by multiplying by 1.678^ 

Con cerning the three galaxies included in the lMancini et al.l 
(|201 0) sample, we note that the r^ values reported in their Ta- 
ble 1 and Figure 4 do not refer to the circularized effective 
radii, but still to the effective semi-major axes a^. Like for the 
other PEGs, here we use properly circularized effective radii 
for objects 307881, 254025, and 217431, using hja = 0.58, 
0.65, and 0.50, respectively (Mancini et al. Erratum in prepa- 
ration). We conclude that the radius of 307881 is in good 
agreement (w ithin 10% and la) with the value reported in 
iMancini et aP (12010). Instead, the effective radii of 254025, 
and 2174 31 are ^ 30%, and ^ 4 5% smaller, respectively, than 
derived inl Mancini et al.l (120101) . though still within ~ 1 .5cr of 
the previous value given the large errors on r^. The Sersic 
indices are all in agreement within ~ 20% (and within ~ 1- 
1 .5a), and the other free parameters of the fit (i.e., magnitude, 
position angle, and axial ratio) are fully consistent (within 2- 
6%, and la) within the published values. All these differences 
can be traced back to the different pixel size (hence PSF sam- 
pling) of the new data and to the different methods of estimat- 
ing the sky background. 

7.2. Stellar mass-size relation 



Figure [2T] shows the resulting effective radii of the pro- 
gram galaxies as a function of the st ellar mass. Fo r com- 
parison, the local relation derived by iNewman et al.l (120121) 
is also shown, together with the subset of z > 1.4 PEGs from 
the compil ation of 465 PEGs with spectroscopic redshift by 
[Damjanov et alj (1201 Ih . For both data sets we have converted 
the stellar masses into those appropriate for the Chabrier IMF 
by usin g correction factors shown in Table 2 of Bemardi et al.l 
(120101) . The majority of these galaxies appear to be more 
compact compared to local PEGs of the same mass, although 
about a third of them lies on, or close to, the local relation 
(within ~ 2a), whereas 7 of them (i.e., ~ 40%) are more than 
3a from the local relation. When only the brighter objects 
with / < 24 are considered, then ~ 50% are classified as com- 
pact. 

Note that the sizes of high-^ galaxies, in particular the 
fainter ones, may have been systematically underestimated 
as shown by simulations (Mancini et al. 2010), an effect that 
may not be confined only to high redshift galaxies. For exam- 
ple, the measured effective radius of M87 can vary by up to a 
facto r of two, depend ing on the depth of the adopted photom- 
etry ("Kormendy et al. 2009). These uncertainties may well 
affect the derived fractions of normal and compact objects in 
our sample. 

Figure[22] shows the effective radius of the 18 galaxies (nor- 
malized to the local relation) vs. their spectr oscopic redshift, 
along with the objects in the compilation of iDamjanov et al.l 
(I2011h . About 40% of our sample show re/re,^=o < 0.3 
while previous studies have found a slightly stronger evolu- 
tion with > 50% of PEGs at z ^ 1 .6 with M > IO^^M^t^ having 
re/re.,=o < 0.3 (e g.,|T ruiillo et al. 2007; Bu itrago et al.ll2008l: 
'Damianov et al.l 120 1 L: iNewman et al.. .2012h . ICooperetaD 
(2012) have shown that PEGs at 0.4 < z < 1.2 in high den- 
sity regions tend to be bigger by up to ~ 25% than those in 
low density regions. If this trend were to continue towards 
higher redshifts then the present PEG sample may be slightly 
biased towards larger sizes compared to those in the general 
field, as our galaxies are preferentially located in overdense 
regions (see Section 4). 

7.3. Relation between Size and Age Indicators 

One effect contributing to the evolution of the mass- 
size relation of PEGs comes from the progressive quench- 
ing of larger sta r-forming gal axies (Ivan der Wei et al.l l2009i 
IValentinuzzi et a l. 2010; Cassa ta et al.ll201 T|). In this scenario, 
one expect that younger, i.e., recently quenched, PEGs have 
size close to the local relation, while older ones, i.e., quenched 
at an earlier time, are more compact. If so, one would ex- 
pect some correlation between the SED-derived age and a de- 
parture from the local mass-size relation. Based on the SED 
ages of 62 spectroscopically identified PEGs at 0.9 < z < 2, 
Saracco et al. ( 201 1) claimed that the normal size high-z ellip- 
tical galaxies have relatively younger stellar population ages 
compared to the compact ones that would show a wider range 
of formation redshifts with a large fraction of them having 
formed at z > 5. In contrast, Whitaker et al.l (1201 2) found that 
younger PEGs at 1.5 < z < 2.0 are more centrally concen- 
trated and may have slightly smaller sizes. To check for this 
effect in our sample, we have compared the size and age indi- 
cators, namely ages from broad-band SED fitting and Dn4000 
indices. 

The left panel of Figure [23] compares the size normalized 
to the local mass-size relation with the age from the SED 
fitting. Besides the present pBzK sample, the SED ages of 
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Figure 21. The stellar mass-size relation of our spectroscopically identi- 
fied pBzKs at 1.4 < z < 1.8 with /gu < 24 (red circles) and /gu > 24 (red 
triangles), where igu is the HST/ACS F814W magnitude recovered by GAL- 
FIT. The orange squares and green diamonds represent passively evolving 
galaxies from the compilation of Damianov et al. ( 2011) with spectroscopic 
redshifts respectively in the range 1.4 <z< 1.8 (the same as in our study) and 
at z > 1.8. The solid line refers to the mean mass-size relation for z = 0.06 
passive galaxies in SDSS from Newman et al. ( 2012) and the gray shaded 
bands show the Icr, 2cr and 3a deviations from such relation. 
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Figure 22. The effective radii of the pBzKs presented here (circles) and 
those of the passive galaxies with spectroscopic redshift from Damianov et al. 
^01 1), normalized to the mass-size relation at z = 0.06 from Newma n et al.i 
(|2012) at the same stellar mass. 



spectroscopically id enti fied PEGs at s imilar redshifts from 
iCimatti et aP (I2008h and lSaracco et alJ (1201 1|) are also sho wn 
for comparison. As pointed out by ISaracco et al.l (1201 lb . in 
their data there appears to be a weak trend with compact 
galaxies showing a wider range of SED ages compared to 
normal- size ones that tend to show younger ages. Splitting 
our sample at log(re/ r^^z^o) = -0.48 (3a deviation from the lo- 
cal mass-size relation), the median of log(Age/yr) are 9.15 
and 9.40 for the normal and compact galaxies, respectively, 
i.e., the compact galaxies appear to be ~ 70% (or about 1 
Gyr) older than the normal ones. When only the objects with 
/814 < 24 are considered, the median log(Age/yr) are 9.10 and 
9.35 for normal and compact objects, respectively. Note that 
the uncertainties in the SED ages are typically comparable to 
these differences. 

The right panel of Figure [23] shows the relation between 
the normalized size and Dn4000, which can be taken as a 
proxy for stellar population age. For compa rison, a com- 
pact quenched galaxy at z = 2.2 dKriek et al.l [2009) is also 
shown. Excluding the object with exceptionally large Dn4000 
(~ 2.75), there seems to be essentially no trend between size 
and Dn4000. The median Dn4000 of the normal and compact 
objects are 1.61 and 1.55, respectively, a difference compara- 
ble to the typical la error in Dn4000 measurement. Restrict- 
ing the sample to the brightest objects, the median Dn4000 of 
normal and compact objects are 1.42 and 1.50, respectively. 

We conclude that in our sample there is no strong correla- 
tion between the size of high redshift PEGs and stellar popu- 
lation age. A Similar conclusion is also found for the brightest 
objects for which the size measurements are more reliable. To 
draw any firm conclusions it is essential to construct a larger 
sample with robustly measured size and Dn4000. 

7.4. Stellar Velocity Dispersion of the pBzK Galaxy 254025 

For the brightest object in the sample (254025) a la upper 
limit to it s stellar velocity disp ersion (a* < 304 km s~^) was 
placed bv lOnoderaeTaH (12010 ) based on our MOIRCS 2009 
data. By adding the data taken in 2010 we have been able to 
upgrade this upper limit to a true measurement of a^ . We used 
a FWHM of 27 A at the observed wavelength for the instru- 
mental profile and assumed the resolution to be constant in 
linear wavelength scale across the wavelength range of inter- 
est. We have restricted the analysis to the 3550-4450 A rest 
frame wavelength range with S/N = 1 1 .5 in 60 km s~^ spectral 
interval, which encompasses features such as CN, Ca II H&K, 
US, the G band, H7 and Fe 4383. We determined a dispersion 
cr* = 270 ± 105 km s"^ from 300 reahzations of pPXF fits that 
consider not only the effect of random noise, but also explore 
the effect of sky residuals and different input parameters in 
pPXF, as described in Section 13.11 Here for a* we take the 
median of the distribution and the la error is determined as 
(cr*, 84 -cr^, 16) /2 where a*, 84 and a^^ie correspond respectively 
to the 84- and 16-percentiles of the probability distribution. 
We emphasize that the quoted error is very conservative, as a 
formal random error of only 32 km s~Ms indicated by single 
pPXF fits. We quote instead an error of 105 km s~^ derived 
above as a more realistic value that includes an estimate of the 
systematic errors. 

Figure [24] compares the structural and dynamical properties 
of the galaxy 254025 with those of local early type galax- 
ies and with those from the small sample of z > 1 PEGs for 
which a stellar velocity dispersion has been measured. Dy- 
namical masses are derived as Mdyn = Sr^al/G where G is 
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Figure 23. Left: The effective radius of passively evolving galaxies with spectroscopic redshifts at z > 1 .4 as a function of their stellar population age from SED 
fitting. Red circles and triangles, green diamonds, and blue squares represent the pBzK sample in our study (depending on their F814W magnitude as in Figure 
1211 . the passively evolving galaxies studied by Saracco et al. ( 2011) at 1.4 < z < 1.8, and those from the GMASS survey (Cimatti et al. 2008) at 1.4 < z < 2. 
Right: Sam e as the left p anel, but expressed as a function of Dn4000. Red symbols are the same as in the left panel, and a yellow square represents a galaxy at 
z = 2.2 from lKriek et al.l ( f200l) . 



the gravitational constant, and we obtained logCMdyn/M©) = 
1 1 .43 ± 0.35 for the galaxy 254025. Figure l24lis an update of 
the similar figure presented in Onodera et al. (2010). We can 
thus confirm that this galaxy closely follows all the local re- 
lations involving the dynamical mass, the effective radius and 
the stellar velocity dispersion within the error bars. Again, 
this is a proof of the coexistence at high redshifts of PEGs 
that appear structurally and dynamically similar to their local 
counterparts with other galaxies that are more compact and 
are correspondingly characterized by a higher velocity dis- 
persion. 

8. NUMBER AND STELLAR MASS DENSITIES 

As our spectroscopic identification of the PEGs with K < 
21 or > 10^^ Mq is nearly complete over the surveyed 
area, we can estimate the abundance of massive PEGs at 
z > 1.4, to within cosmic variance errors. The redshift range 
covered by the present sample 1.43 < z ^ 1.83 over the sur- 
vey area of 66.2 arcmin^ contains a comoving volume of 
^ 82,000 Mpc^ Thus, from the 16 objects with ^ < 21 we 
estimate a comoving number density of 1 .9 x 10~^ Mpc~^, and 
the corresponding stellar mass density is 3 .0 x lO^M© Mpc~^, 
with an uncertainty of ~ 25% from Poisson statistics. These 
should be regarded as upper limits since the field was chosen 
for its overdensity of PEGs. However they are in agreement 
with estimates derived from a complete sample of galaxies in 
the Hubble Ultra Deep Field of 3.4 x lO'^ Mpc"^ and 1.9- 
4.6 X 10'^ M cT) Mpc~^ (after con verting to Chabrier IMF), as 
measured by Dadd i et al.l (|2005|) for spectroscopically identi- 
fied PEGs at 1.39 < z < 2.00 down to the limit K = 23, al- 
though their explo red co moving volume is ~ 4 times smaller. 

iBrammer et aH (|2011l) have derived number and mass den- 
sities of PEGs at 0.4 < Zphot < 2.2 from the NEWFIRM 
Medium Band Survey (NMBS), covering an area of '^0.4 
deg^, and for the mean redshift of our galaxies (z = 1.69) 



found (1.0 ±0.3) X lO'^ Mpc"^ and (1.6 ±0.4) x 10^ M© 
Mpc~^, respectively. Our estimates are somewhat higher, as 
expected given the mentioned overdensity bias. 

Comparing to the corresponding local values, respectively 
-1x1 0"4 Mpc"^ and - 8 X lO^M© Mpc"^ from lBaldrv et al.l 
(12004)) . taken at face value our result suggests that PEGs at 
1 .4 < z < 1 .8 account for — 50% of the number density and — 
35-40% of the stellar mass density of local PEGs with > 



9. SUMMARY AND CONCLUSIONS 



Mapping the 3D space distribution of passively evolving 
galaxies at high redshift is an endeavour of critical importance 
in the broader context of galaxy evolution and observational 
cosmology. However, it is perhaps also the most challenging 
observationally, as it requires many hundreds of hours of 8m- 
class telescope time with current instrumentation. Emission 
lines are absent (or very weak) in the spectra of PEGs, and 
therefore measuring redshifts and velocity dispersions must 
rely uniquely on absorption lines. With spectrometers work- 
ing at optical wavelengths, for such objects at z ^ 1.4 this 
is possible using absorption lines over an intrinsically very 
faint UV continuum, which requires extremely long integra- 
tion times. Using near-IR, multiobject spectrometers offers a 
potentially attractive alternative, as the strongest spectral fea- 
tures over a strong continuum become accessible. However, 
ground based, infrared spectroscopy of faint, high redshift 
objects is affected by a higher sky and thermal background 
which largely limits this advantage. Using the MOIRCS in- 
strument at the Subaru telescope we have started a pilot exper- 
iment aimed at studying a sample of bright PEGs at z ^ 1.4 
while assessing the feasibility of a more ambitious project. 

With this objective, we have first identified the strongest 
concentrations of high-z candidate PEGs over the COSMOS 
field, and focused our observational effort over one such con- 
centration and its immediate vicinity. A sample of 34 objects 
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Figure 24. A comparison of the structural/dynamical properties of high-redshift early-type galaxies with measured velocity dispersions (symbols with error 
bars) with those from SDSS elliptical galaxies at z ^ 0.06 (shown as gray scale), (a): The stellar velocity dispersions vs. the dynamical mass, (b): Comparison 
between the dynamical and the stellar masses; the diagonal dashed line refers to equal mass and the dotted lines show the range around it by a factor of 2. (c): 
The effective radii vs. the dynamical masses. {d)\ The effective radii vs. the stellar velocity dispersions. In all panels, the symbols refer to: the galaxy 254025 
from the present study (red filled circle); an object at z = 2.19 (green filled circle; van Dokkum et al. 2009); two GMASS galaxies at z = 1.41 and the stacked 
properties of the GMAS S galaxies at 1.6 < z < 2.0 (blue squares and yello w triangle, respective ly; Cappellari et al...2009.) : an object at z = 1.80 (cyan circle; 
Ivan de Sande et al.|[201 II) ; and objects at 1.05 < z < 1.60 (small pink circles; iNewman et aT1l201QI) . 



has been observed, with integration times ranging from 280 
to 550 minutes, and the continuum has been detected in 18 of 
them. The redshift of these detected galaxies are concentrated 
in three spikes at z = 1.43, 1.53 and 1.82, and the last two 
galaxies have quite similar redshifts around 1.67. This spec- 
troscopically confirms the strong clustering of high redshift 
PEGs, previously hinted only from th eir 2-point correla tion 
function (e.g., Kong et al. 2006; Mc Cracken et al.ll201Qh . A 
comparison with COSMOS photometric redshifts shows that 
2/3 of them agree closely with the spectroscopic redshifts. 



whereas the other ~ 1 /3 of them are systematically underes- 
timated by up to ~ 25%. Thus, these photometric redshifts 
can lead to a sizable underestimate of the volume density of 
PEGs at z > 1.4. However, we show that photometric redshift 
training with these spectroscopic redshifts can significantly 
improve their performance for high-z PEGs. 

Having fixed the redshifts to their spectroscopic values, we 
have explored the stellar population content of these galaxies, 
measuring their masses, luminosity-weighted ages and metal- 
licities using various SED-fitting techniques. Stellar masses 
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range from ~ 4 x 10^^ to ~ 4 x 10^^ M©, ages from ~ 1 to ~ 3 
Gyr, whereas near-solar metallicity is found for most galaxies. 

The Lick H^p and the Dn4000 indices allowed us to inves- 
tigate further the stellar population content and evolutionary 
stage of the sample galaxies. We find that about 40% of the 
pBzKs show large Dn4000 (around 2) which indicates that 
these galaxies are already quenched long ago compared to 
their timescale of star formation. On the other hand, about 
60% of them have Dn4000 values close to those of young 
star-forming or post-starburst galaxies, which indicates that 
these galaxies have been quenched more recently, given the 
lack of the star-forming activity (except for one object with a 
mid-IR detection). This is also confirmed for the four most 
massive and brightest objects and the composite spectrum 
with measurable absorption line. Two of them (as well as 
the composite spectrum, see below) show moderate Dn4000 
and H^F, consistent with the early stages of a post-starburst 
galaxy where star formation has been quenched very recently. 
The other two bright galaxies show both large Dn4000 and 
H^F and appear to be already passive, possibly with a recent 
episode of star-formation, though the combination of two in- 
dices seems rather extreme. 

The composite spectrum of the 17 galaxies exhibits sev- 
eral well-detected spectral features, notably CN, Ca II H&K, 
the G band, Mgb and various iron lines. A fit to 

this composite spectrum indicates that it is dominated by old, 
metal-rich stellar populations with a mass-weighted age of 
~ 3.5 Gyr. The fit is somewhat dependent on the allowed 
age and metalhcity ranges, but the existence of stellar popu- 
lation younger than ~ 400 Myr is always excluded by the fits, 
which is confirmed by the use of another fitting code adopting 
a different input SSP library. 

The structural parameters, (circularized) effective radius 
and Sersic index, have been measured for these galaxies over 
the HST/ ACS /-band image of the COSMOS field. A major- 
ity of the 1 8 galaxies have effective radii somewhat lower than 
equally massive galaxies in the local Universe, with some up 
to three times lower. However, almost about a third of the 
galaxies appear to be on, or close to, the local r^-M^ rela- 
tion, confirming that compact and apparently normal PEGs 
coexist at these redshifts. 

Finally, the stellar velocity dispersion of the brightest 
galaxy in the sample was measured as a* = 270 ±105 km s~^ . 
Combined with its effective radius and stellar mass, this value 
of cr* confirms that this z = 1 .82 galaxy closely follows all the 
local structural and dynamical relations for PEGs. 

The nearly complete identification of PEGs at the top end 
of their stellar mass function at the observed redshift range 
allows us to estimate their number and stellar mass densities, 
still with the proviso that they refer to an overdensity region. 
Thus, the PEGs with M* > 10^^ M© at 1 .4 > z > 1 .8 we found 
in the studied field correspond to about 50% of the number 
and about 35-40% of the stellar mass densities of local PEGs 
with similar stellar mass. 

Without counting the time lost to technical problems and 
bad weather, this project has used 4-5 nights of Subaru 
telescope time. Over the whole 2 square degree COS- 
MOS field one can roughly count ~ 20 concentrations of 
candidate high-z PEGs similar to the one explored in the 
present study. Thus, some 80-100 nights would be re- 
quired to map all such concentrations in a similar fashion 
with Subaru/MOIRCS. Obviously, even more would be re- 
quired to map the distribution of all massive PEGs over the 



entire COSMOS field. A coordinated effort with the up- 
coming near-IR multi-object spectrographs at 8-10 m class 
telescopes (such as Gemini/FLAMINGOS-II, VLT/KMOS, 
and Keck/MOSFIRE) would be required to make such large 
project feasible in the near future. 
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